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43 The Laplace Transform: Basic Definitions 
and Results 

Laplace transform is yet another operational tool for solving constant coeffi¬ 
cients linear differential equations. The process of solution consists of three 
main steps: 

• The given “hard” problem is transformed into a “simple” equation. 

• This simple equation is solved by purely algebraic manipulations. 

• The solution of the simple equation is transformed back to obtain the so¬ 
lution of the given problem. 

In this way the Laplace transformation reduces the problem of solving a dif¬ 
ferential equation to an algebraic problem. The third step is made easier by 
tables, whose role is similar to that of integral tables in integration. 

The above procedure can be summarized by Figure 43.1 



Figure 43.1 


In this section we introduce the concept of Laplace transform and discuss 
some of its properties. 

The Laplace transform is defined in the following way. Let f(t ) be defined 
for t > 0. Then the Laplace transform of /, which is denoted by £[/(£)] 
or by F(s), is defined by the following equation 

pT poo 


f(t)e st dt = / /(£)e st dt 


£[/(*)] = F ( s ) = l im 

T—► oc 

The integral which defined a Laplace transform is an improper integral. An 
improper integral may converge or diverge, depending on the integrand. 
When the improper integral in convergent then we say that the function /(£) 
possesses a Laplace transform. So what types of functions possess Laplace 
transforms, that is, what type of functions guarantees a convergent improper 
integral. 


Example 43.1 

Find the Laplace transform, if it exists, of each of the following functions 


(a) f(t) = e at ( b ) f(t) = 1 (c) f(t) = t ( d ) f(t) = e * 2 
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Solution. 

(a) Using the definition of Laplace transform we see that 



But 




T 

l_ e -(s-°)T 


if s = a 
if s ^ a. 


For the improper integral to converge we need s > a. In this case, 


C[e at ] = F(s) 


1 

-, s > a. 

s — a 


(b) In a similar way to what was done in part (a), we find 


poo pi 

£[1] = / e~ st dt = lim 

Jo o 


e dt = -, s > 0. 
s 


(c) We have 


C[t] — / te st dt = 


te st e st 


= s>0 - 


(d) Again using the definition of Laplace transform we find 


poo 

C[e t2 } = / e t2 ~ st dt 

Jo 

If s < 0 then t 2 — st > 0 so that > 1 and this implies that / 0 °° e t2 ~ st dt > 
J 0 °° . Since the integral on the right is divergent, by the comparison theorem 
of improper integrals (see Theorem 43.1 below) the integral on the left is also 
divergent. Now, if s > 0 then J 0 °° e^-^dt > J s °° dt. By the same reasoning 
the integral on the left is divergent. This shows that the function f(t) = e* 2 
does not possess a Laplace transform ■ 


The above example raises the question of what class or classes of functions 
possess a Laplace transform. Looking closely at Example 43.1(a), we notice 
that for s > a the integral J 0 °° e~^ s ~ a2>t dt is convergent and a critical compo¬ 
nent for this convergence is the type of the function f(t). To be more specific, 
if f{t) is a continuous function such that 

\f(t)\<Me at , t > C (1) 
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where M > 0 and a and C are constants, then this condition yields 


rC 


f(i)e~ st dt < / f(t)e~ st dt + M / e^ s ~ a)t dt. 


ic 


Since /(f) is continuous in 0 < t < C, by letting A = max{|/(f)| : 0 < t < C} 
we have 


J f(t)e st dt < A J e st dt = A ^-^ < oo. 

On the other hand, Now, by Example 43.1(a), the integral J^° e^ (s_a ' ) *dt is 
convergent for s > a. By the comparison theorem of improper integrals (see 
Theorem 43.1 below) the integral on the left is also convergent. That is, /(f) 
possesses a Laplace transform. 

We call a function that satisfies condition (1) a function with an exponential 
order at infinity. Graphically, this means that the graph of /(f) is contained 
in the region bounded by the graphs of y — Me at and y = —Me at for f > C. 
Note also that this type of functions controls the negative exponential in the 
transform integral so that to keep the integral from blowing up. If C = 0 
then we say that the function is exponentially bounded. 

Example 43.2 

Show that any bounded function /(f) for f > 0 is exponentially bounded. 

Solution. 

Since /(f) is bounded for f > 0, there is a positive constant M such that 
|/(f)| < M for all f > 0. But this is the same as (1) with a = 0 and C — 0. 
Thus, fit ) has is exponentially bounded ■ 

Another question that comes to mind is whether it is possible to relax the 
condition of continuity on the function /(f). Let’s look at the following situ¬ 
ation. 

Example 43.3 

Show that the square wave function whose graph is given in Figure 43.2 
possesses a Laplace transform. 
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Figure 43.2 

Note that the function is periodic of period 2. 

Solution. 

Since f{t)e~ st < e~ st , we have / 0 °° f(t)e~ st dt < J 0 °° e~ st dt. But the integral on 
the right is convergent for s > 0 so that the integral on the left is convergent 
as well. That is, £[/(£)] exists for s > 0 ■ 

The function of the above example belongs to a class of functions that we 
define next. A function is called piecewise continuous on an interval if 
the interval can be broken into a finite number of subintervals on which the 
function is continuous on each open subinterval (i.e. the subinterval without 
its endpoints) and has a finite limit at the endpoints (jump discontinuities 
and no vertical asymptotes) of each subinterval. Below is a sketch of a 
piecewise continuous function. 



Figure 43.3 

Note that a piecewise continuous function is a function that has a finite 
number of breaks in it and doesnt blow up to infinity anywhere. A function 
defined for t > 0 is said to be piecewise continuous on the infinite in¬ 
terval if it is piecewise continuous on 0 < t < T for all T > 0. 


Example 43.4 

Show that the following functions are piecewise continuous and of exponential 
order at infinity for t > 0 
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(a) f{t) = t n (b) fit) = t n sin at 


Solution. 

(a) Since e t = YJ™=o we l iave f n < rile*. Hence, t n is piecewise con¬ 

tinuous and exponentially bounded. 

(b) Since |f n sinaf| < n!e/ we have t n sin at is piecewise continuous and ex¬ 
ponentially bounded ■ 

Next, we would like to establish the existence of the Laplace transform for 
all functions that are piecewise continuous and have exponential order at 
infinity. For that purpose we need the following comparison theorem from 
calculus. 

Theorem 43.1 

Suppose that fit) and gift) are both integrable functions for all t > t 0 such 
that |/(f) | < | git) for f > f 0 . If f^ g(t)dt is convergent, then f(t)dt is 
also convergent. If, on the other hand, f(t)dt is divergent then f™ f(t)dt 
is also divergent. 

Theorem 43.2 (Existence) 

Suppose that fit) is piecewise continuous on t > 0 and has an exponential 
order at infinity with |/(t)| < Me at for t > C. Then the Laplace transform 

poo 

F(s) = / f(t)e~ st dt 
Jo 

exists as long as s > a. Note that the two conditions above are sufficient, but 
not necessary, for F(s) to exist. 

Proof. 

The integral in the definition of F(s) can be splitted into two integrals as 
follows 

poo pc poo 

/ f(t)e~ st dt= / f(t)e~ st dt + / f(t)e~ st dt. 

Jo Jo Jc 

Since fit) is piecewise continuous in 0 < t < C, it is bounded there. By 
letting A = max{|/(f)| : 0 < f < C} we have 

f c f\ e~ sC \ 

f(t)e st dt < A J e st dt = A (- J < oo. 
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Now, by Example 43.1(a), the integral f ( f{t)e~ st dt is convergent for s > a. 
By Theorem 43.1 the integral on the left is also convergent. That is, f(t) 
possesses a Laplace transform i 

In what follows, we will denote the class of all piecewise continuous func¬ 
tions with exponential order at infinity by VS. The next theorem shows that 
any linear combination of functions in VS is also in VS. The same is true for 
the product of two functions in VS. 

Theorem 43.3 

Suppose that f(t ) and g(t ) are two elements of VS with 

\f(t)\ < Mie ait , t > Ci and \g(t)\ < M 2 e ait , t > C 2 . 

(i) For any constants a and (3 the function af(t) + (3g(t) is also a member of 
VS. Moreover 

C[af{t) + pg(t)\ = aC[f(t)} + pC[g{t)\. 

(ii) The function h(t) = f(t)g(t) is an element of VS. 

Proof. 

(i) It is easy to see that af(t ) + f3g(t) is a piecewise continuous function. 
Now, let C = C\ + C 2 , a = max{ai, a 2 }, and M = |cr|Afi + \/3\M 2 . Then for 
t > C we have 

\<xf(t)+Pg(t)\ < \a\\f(t)\ + \P\\g(t)\ < \a\M ie a ^ + \/3\M 2 e a ^ < Me at . 

This shows that a.fit) + /3g(t ) is of exponential order at infinity. On the 
other hand, 

C[af(t) + Pg(t)] = limT^oo J^[af{t) + fig(t)]dt 

= a liniT^oo f 0 T f {t)dt + P lim^oo / Q T g(t)dt 
= aC[f(t)]+pC[g(t)] 

(ii) It is clear that h(t) = f(t)g(t) is a piecewise continuous function. Now, 
letting C = Ci + C 2 , M = M l M 2l and a = ai + a 2 then we see that for t > C 
we have 


\h(t)\ = |/(t)|| 0 (t)| < MiM 2 e {ai+a2)t = Me at . 



Hence, h(t) is of exponential order at infinity. By Theorem 43.2 , C[h(t)] 
exists for s > a M 

We next discuss the problem of how to determine the function /(f) if F(s) 
is given. That is, how do we invert the transform. The following result on 
uniqueness provides a possible answer. This result establishes a one-to-one 
correspondence between the set V£ and its Laplace transforms. Alterna¬ 
tively, the following theorem asserts that the Laplace transform of a member 
in VS is unique. 

Theorem 43.4 

Let /(f) and g(t) be two elements in VS with Laplace transforms F(s) and 
G(s) such that F(s) = G(s) for some s > a. Then /(f) = g(t) for all t > 0 
where both functions are continuous. 

The standard techniques used to prove this theorem( i.e., complex analysis, 
residue computations, and/or Fourier’s integral inversion theorem) are gen¬ 
erally beyond the scope of an introductory differential equations course. The 
interested reader can find a proof in the book ” Operational Mathematics” 
by Ruel Vance Churchill or in D.V. Widder ’’The Laplace Transform”. 

With the above theorem, we can now officially define the inverse Laplace 
transform as follows: For a piecewise continuous function / of exponential 
order at infinity whose Laplace transform is F. we call / the inverse Laplace 
transform of F and write / = £ _ 1 [.F(s)]. Symbolically 

/(*)= ■C-'If’M] <=*• F(s) = £[/(*)]. 


Example 43.5 

Find £ _1 (^j) , s > 1. 

Solution. 

From Example 43.1(a), we have that C[e at ] = V—, s > a. In particular, for 
a = 1 we find that £[e*] = s > 1. Hence, C~ l (y^y) = e*, t > 0 |. 

The above theorem states that if fit) is continuous and has a Laplace trans¬ 
form F(s), then there is no other function that has the same Laplace trans¬ 
form. To find £ _1 [E(s)], we can inspect tables of Laplace transforms of 
known functions to find a particular /(f) that yields the given F(s). 

When the function /(f) is not continuous, the uniqueness of the inverse 
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Laplace transform is not assured. The following example addresses the 
uniqueness issue. 

Example 43.6 

Consider the two functions /(f) = h(t)h (3 — f) and g(t) = h(t ) — h(t — 3). 

(a) Are the two functions identical? 

(b) Show that £[/(£)] = C[g(t). 

Solution. 

(a) We have 


and 


So the two functions 
(b) We have 

f 3 1 — e~ 3s 

£[/(*)] = C\g(t)] = / e st dt = -, s > 0. 

Jo s 

Thus, both functions /(f) and g[t) have the same Laplace transform even 
though they are not identical. However, they are equal on the interval(s) 
where they are both continuous ■ 

The inverse Laplace transform possesses a linear property as indicated in 
the following result. 

Theorem 43.5 

Given two Laplace transforms F(s) and G(s) then 

£- ! K(s) + bG(s)] = a£ _ 1 [F(s)] + 6 £- 1 [G(s)] 
for any constants a and b. 

Proof. 

Suppose that £[/(£)] = F(s) and C[g{t)] = G(s). Since C[af(t) + bg(t)] = 
aC[f(t)\ + bC[g{t)] = aF(s ) + bG(s) we have £ _ 1 [oF(s) + bG(s)] = a fit) + 
bg(t) = a£ _ 1 [F(s)] + 6 £ _ 1 [G(s)] ■ 


m = 


1, 0 < t < 3 
0, f > 3 


g(t) = 


i, 

o, 


0 < t < 3 
t > 3 


are equal for all t 3 and so they are not identical. 
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Practice Problems 


Problem 43.1 

Determine whether the integral J 0 °° j^dt converges. If the integral con¬ 
verges, give its value. 

Problem 43.2 

Determine whether the integral J 0 °° j^pdt converges. If the integral con¬ 
verges, give its value. 

Problem 43.3 

Determine whether the integral J 0 °° e -< cos (e _t )df converges. If the integral 
converges, give its value. 

Problem 43.4 

Using the definition, find C[e 3t ], if it exists. If the Laplace transform exists 
then find the domain of F(s). 

Problem 43.5 

Using the definition, find C[t — 5], if it exists. If the Laplace transform exists 
then find the domain of F(s). 

Problem 43.6 

Using the definition, find ], if it exists. If the Laplace transform 

exists then find the domain of F(s). 

Problem 43.7 

Using the definition, find C[(t — 2) 2 ], if it exists. If the Laplace transform 
exists then find the domain of F(s). 

Problem 43.8 

Using the definition, find £[/(£)], 
then find the domain of F(s). 

m = { 


if it exists. If the Laplace transform exists 

0, 0 < t < 1 

t - 1, t > 1 
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Problem 43.9 

Using the definition, find £[/(£)], if it exists. If the Laplace transform exists 
then find the domain of F(s). 

( 0, 0 < t < 1 

/(f) = l f-1, 1 < t < 2 
[ 0, t > 2. 

Problem 43.10 

Let n be a positive integer. Using integration by parts establish the reduction 
formula 

j t n e~ st dt = + 1 ! J t n -\- st dt, s > 0. 

Problem 43.11 

For s > 0 and n a positive integer evaluate the limits 
lirrp^o t n e~ st (b) linp^oo t n e~ st 


Problem 43.12 

(a) Use the previous two problems to derive the reduction formula for the 
Laplace transform of /(f) = f n , 

C[t n ] = '^C{t n -\ s > 0. 

(b) Calculate C[t k ], for k = 1,2, 3,4, 5. 

(c) Formulate a conjecture as to the Laplace transform of /(f), t n with n a 
positive integer. 

From a table of integrals, 

f e au sin /3udu = 
f e au cos fdudu = 

Problem 43.13 

Use the above integrals to find the Laplace transform of fit) = costat, if it 
exists. If the Laplace transform exists, give the domain of F(s). 

Problem 43.14 

Use the above integrals to find the Laplace transform of /(t) = sin tat, if it 
exists. If the Laplace transform exists, give the domain of F(s ). 


au & sin (3u—(3 sin (3u 
e ^£2 

au ol cos /3n+/3 sin flu 
e a 2 +(3 2 
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Problem 43.15 

Use the above integrals to find the Laplace transform of /(f) = cos cu(t — 2), 
if it exists. If the Laplace transform exists, give the domain of F(s). 

Problem 43.16 

Lise the above integrals to find the Laplace transform of /(f) = e 3t sin f, if it 
exists. If the Laplace transform exists, give the domain of F(s). 

Problem 43.17 

Use the linearity property of Laplace transform to find C[5e~ n + t + 2e 2t ]. 
Find the domain of F(s). 

Problem 43.18 

Consider the function /(f) = tanf. 

(a) Is /(f) continuous on 0 < f < oo, discontinuous but piecewise contin¬ 
uous on 0 < t < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 

Problem 43.19 

Consider the function /(f) = f 2 e _t . 

(a) Is /(f) continuous on 0 < f < oo, discontinuous but piecewise contin¬ 
uous on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 

Problem 43.20 

t 2 

Consider the function /(f) = J t+1 . 

(a) Is /(f) continuous on 0 < f < oo, discontinuous but piecewise contin¬ 
uous on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 

Problem 43.21 

Consider the floor function /(f) = |/J, where for any integer n we have 
[_fj = n for all n < t < n + 1. 

(a) Is f(t ) continuous on 0 < f < oo, discontinuous but piecewise contin¬ 
uous on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 
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Problem 43.22 

Find (£) . 

Problem 43.23 

Find £-'(-?+i?l)- 

Problem 43.24 

Find C~ l (^2 + ■ 
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44 Further Studies of Laplace Transform 

Properties of the Laplace transform enable us to find Laplace transforms 
without having to compute them directly from the definition. In this sec¬ 
tion, we establish properties of Laplace transform that will be useful for 
solving ODEs. 


Laplace Transform of the Heaviside Step Function 

The Heaviside step function is a piecewise continuous function defined by 


h{t) = 


1, t > 0 
0, t < 0 


Figure 44.1 displays the graph of h(t). 


; 

y 

l 




t 


Figure 44.1 


Taking the Laplace transform of h{t) we find 

/»oo 

C[h(t)] = / h(t)e~ st dt = 


/»oo 

r n — st~\ 

/ e~ st dt = 

C 

Jo 

s 


= -, s > 0. 


A Heaviside function at a > 0 is the shifted function h(t — a) (a units to the 
right). For this function, the Laplace transform is 


roo 

r e~ st ~\ 

/ e~ st dt = 

e 

la 

s 


C[h{t — a)] = / h(t — a)e st dt = 


Laplace Tranform of e at 

The Laplace transform for the function f(t) = e at is 


-, s > 0. 


C[e at ] = / e- {s ~ a)t dt = 


e ~(s-a)t 


s — a 


s — a 


-, s > a. 
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Laplace Tranforms of sin at and cos at 

Using integration by parts twice we find 

£[sin at] = J 0 °° e~ st sin atdt 


e st sin at ae st cos at 


^ r e st s ^ n a ^ 


— ^£ [sin at] 


s 2 +a 2 


j £ [sin at] = 


a 

<5 2 


£ [sin at] = 

A similar argument shows that 

£ [cos at] = 


2 -L 2 , s > 0 


s z +a 


s 2 + a 2 


, s > 0. 


Laplace Transforms of cosh at and sinh at 

Using the linear property of £ we can write 

£[cosh at] = \ (£[e at ] + £[e -Qi ]) 

= b (— + tr), s > |, 

2 \s—a s+a/ ’ 1 


A similar argument shows that 

£ [sin at] = 


9 9 ’ 

— cr 


s > a 


s > a . 


Laplace Transform of a Polynomial 

Let n be a positive integer. Using integration by parts we can write 


t n e~ st dt = 


t n e~ st 


n 
H— 
s 


t n ~ 1 e~ st dt. 


By repeated use of L’Hopital’s rule we find Hindoo t n e st = lim^oo 
for s > 0. Thus, 

£[t n ] = -^C[t n ~\ s > 0. 


= 0 
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Using induction on n — 0,1, 2, • • • one can easily eastablish that 


C[t n ] 


n\ 

s n+1 ’ 


s > 0. 


Using the above result together with the linearity property of C one can find 
the Laplace transform of any polynomial. 

The next two results are referred to as the first and second shift theorems. 
As with the linearity property, the shift theorems increase the number of 
functions for which we can easily find Laplace transforms. 


Theorem 44.1 (First Shifting Theorem) 

If f{t) is a piecewise continuous function for t > 0 and has exponential order 
at infinity with \f{t)\<Me at , t > C , then for any real number a we have 

C[e at f(t)] = F(s — a), s > a + a 

where £[/(£)] = F(s ). 

Proof. 

From the definition of the Laplace transform we have 

poo poo 

C[e at f(t )] = / e~ st e at f(t)dt = / e~^- a)t f{t)dt. 

Jo Jo 

Using the change of variable (3 = s — a the previous equation reduces to 

poo poo 

C[e at f(t)}= / e~ st e at f{f)dt= / e~ pt f(t)dt = F(fi) = F(s-a), s > a+a m 
Jo Jo 

Theorem 44.2 (Second Shifting Theorem) 

If f{t) is a piecewise continuous function for t > 0 and has exponential order 
at infinity with \f(t)\ < Me at , t > C, then for any real number a > 0 we 
have 

C[f(t — a)h(t — a)] = e~ as F(s), s > a 
where C[f{t)] = F(s) and h(t) is the Heaviside step function. 

Proof. 

From the definition of the Laplace transform we have 

poo poo 

C[f(t — a)h(t — a)] = / f(t — a)h(s — a)e~ st dt = / f(t — a)e~ st dt. 

Jo Ja 
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Using the change of variable (3 = t — a the previous equation reduces to 
C[f{t - a)h(t - a)] = / 0 °° f(/3)e~ s ^ +a) d/3 

— e~ sa f(/3)e~ sl5 df3 = e~ sa F(s), s > a M 

Example 44.1 

Find 


(a) C[e 2t t 2 ] (b) £[e 3 *cos2f] (c) C 1 [e 2t s 2 


Solution. 

(a) By Theorem 44.1, we have C[e 2t t 2 ] = F(s — 2) where C[t 2 } — jl — 

F(s), s > 0. Thus, C[e 2t t 2 } = s > 2. 

(b) As in part (a), we have £[e 3t cos2f] = F(s — 3) where £[cos2f] = F(s — 3). 
But £[cos2f] = ppj, s > 0. Thus, 

£[e 3t cos2f] = (g ^ 3)2 3 +4 , s>3 


(c) Since C[t] = by Theorem 44.2, we have 

-2t 

- = C[(t-2)h(t-2)}. 
s z 

Therefore, 


C - 1 



(■t - 2 )h{t - 2) 


0, 0 < t < 2 

t- 2, t >2 u 


The following result relates the Laplace transform of derivatives and integrals 
to the Laplace transform of the function itself. 


Theorem 44.3 

Suppose that /(f) is continuous for t > 0 and /'(f) is piecewise continuous 
of exponential order at infinity with |/'(f)j < Me at , t > C Then 


(a) /(f) is of exponential order at infinity. 

(b) £[/'(f)] = s£[/(f)l — /(0) = sF(s ) — /(0), s > max{a,0} + 1. 

(c) £[/"(f)] = s 2 C[f(t)} - sf{ 0) - /'(0) = , 2 F( S ) - s/(0) - /(0), a > 
max{a, 0} + 1. 


( d ) £ fo f ( u )du 




s > max{a, 0} + 1. 
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Proof. 

(a) By the Fundamental Theorem of Calculus we have f(t) = /(0)~ f'{u)du. 
Also, since f is piecewise continuous, < T for some T > 0 and all 

0 < t < C. Thus, 


i/eoi 


/(°) - Jo f'( u )du 


l/(°) - Jo f( u ) du ~ Sc f(u)du 


< 

Note that if a > 0 then 


|/(0) | + TC + M fS e au du 


and so 


If a = 0 then 


and therefore 




a 


e at 

< — 
a 


1/(01 < [|/( 0 )| + rc + ^]e 

Lb 


at 


e au du = t-C 


ic 


\m\ < |/(0)| +TC + M(t -C)< (|/(0)| +TC + M)e t . 
Now, if a < 0 then 


r e au du= ~(e at — e aC ) < — 


'C 


a 


a 


so that 


It follows that 


1/(01 <(l/(0)| + rc + ^)e 


\f(t)\<Ne u , t> 0 


where b = max{a, 0} + 1. 

(b) From the definition of Laplace transform we can write 


£[/'(*)] = lim / f(t)e st dt. 

A -^OO /n 
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Since f'(t) may have jump discontinuities at An in the interval 

0 < t < A, we can write 


nvt 


Integrating each term on the RHS by parts and using the continuity of f(t) 
to obtain 


So 1 f'(t)e st dt = f(ti)e stl - /(0) + s/J 1 f(t)e st dt 

f{t)e~ st dt = f(t 2 )e ~ st2 - / (ti)e~ stl + s f(t)e~ st dt 


ItH-ifi *) 6 Stdt = f(tN)e stN -f(t N - i)e stN ^ f{t)e st dt 


tN 


ftt /'(*) e Stdt = /(^) e aA “ f(t N )e stN + s f A f(t)e st dt 

Also, by the continuity of f(t) we can write 


f[t)e st dt= f f{t)e st dt+ f f(t)e st dt + ---+ f f(t)e st dt. 

JO J t\ J tjy 


Hence, 


f(t) e ~ st dt = f(A)e~ sA - /(0) + s / f(t)e~ st dt. 


Since f(t) has exponential order at infinity,lim^oo f(A)e sA = 0. Hence, 

C[f(t)} = sC[f(t)}-f( 0). 


(c) Using part (b) we find 

c[nt)} = sc[f(t)} - n o) 

s{sF{s) - /(0)) - /'(0) 

= s 2 F(s ) — sf( 0) — /'(0), s > max{a, 0} + 1 


(d) Since A ( f(u)du j = f(t), by part (b) we have 


F(s) = C[f{t )] = sC \ / f{u)du 
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and therefore 


C I f(u)du = — — — = — — , s > max{a, 0} + 1 ■ 

.Jo \ s s 

The argument establishing part (b) of the previous theorem can be extended 
to higher order derivatives. 

Theorem 44.4 

Let /(£), ■ ■ ■ , be continuous and f^ n \t) be piecewise continu¬ 

ous of exponential order at infinity with \f^ n \t)\ < Me at , t > C. Then 

£[/«”>(«)] = S ”£[/(t)]- S "-'/(0)- S "- 2 /'(0)-/<”-'>(0), s > max{a, 0}+l. 

We next illustrate the use of the previous theorem in solving initial value 
problems. 

Example 44.2 

Solve the initial value problem 

y"-4y' + 9y = t, y(0) = 0* j/(0) = 1. 

Solution. 

We apply Theorem 44.4 that gives the Laplace transform of a derivative. By 
the linearity property of the Laplace transform we can write 

C[y"} — 4£[j/'] + 9£[y\ = C[t]. 

Now since 

C[y"} = s 2 C[y} - sy{0) - y'(0) = s 2 Y(s) - 1 
C[y'} = sY(s) - y(0) = sY(s) 

m = $ 

where C[y) = E(s), we obtain 

s 2 l"(s) - 1 - 4sF(s) + 9T(s) = 4. 

s z 

Rearranging gives 

(s 2 -As + 9)Y(s) = 

s z 
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Thus, 


and 


Y{s) = 

s 2 (s 2 

y{t) = C~ l 

_s 2 (s 2 


s 2 + 1 


- 4s + 9) 
? 2 + 1 


In the next section we will discuss a method for finding the inverse Laplace 
transform of the above expression. 


Example 44.3 

Consider the mass-spring oscillator without friction: y" + y = 0. Suppose 
we add a force which corresponds to a push (to the left) of the mass as it 
oscillates. We will suppose the push is described by the function 

f(t ) = —h(t — 2n) + u(t — (27t + a)) 

for some a > 2tt which we are allowed to vary. (A small a will correspond 
to a short duration push and a large a to a long duration push.) We are 
interested in solving the initial value problem 

y" + y = f(t), y ( o) = l, j/'(o) = o. 


Solution. 

To begin, determine the Laplace transform of both sides of the DE: 

£[y" + y\ = £[/(*)] 


or 


Thus, 


Now since 


s 2 Y - WO) - WO) + Y(s) = —~e~ 27rs + -e _(27r+a)s . 

s s 


e —(2 7r+a)s g—: 2 tts s 

Y (s) = / n . . X-/ o . . X + 


1 _ 1 


s(5 2 + 1) s s 2 + 1 


s(s 2 + 1) s(s 2 + 1) s 2 + 1 

we see that 


(2n-\-a)s 

1 s 

—2-its 

ri . i 


Co | 

CO 

to 

+ 

h- 1 

o 

CO | 

CO 

to 

+ 

I— 1 

t_ 


+ 


s 2 + 1 
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and therefore 


y(t) = h{t- (27T + a)) [C 1 (i - ^)] (t- ( 2 7T + a)) 

h(t - 2tt) [£ _1 Q - -^)] (t - 2 tt) + cost 
= hit — (27t + a))[1 — cos (t — (27 r + a))] — u(t — 2n) [1 — cos (t — 27t)] 
+ cos t ■ 


We conclude this section with the following table of Laplace transform pairs. 


m 

h(t ) = 


1, t > 0 
0, t< 0 


f\ n = 1, 2, • • • 


sin {out) 
cos (c ot) 
sinh (c ot) 


F(s) 

s 

O 

A 

s ’ 


n\ 

s n+l •) 

s > 0 

S 

s > a 

s—a ’ 


UJ 

S 2 +Cd 2 

, s > 0 

S 

5 2 +CJ 2 

, s > 0 

o~ 2 , s > cu 

S z —U) z 7 1 1 


cosh (t ot) 


S 

s 2 —uj 


2 5 


S > \UJ 


e at f{t ), with \f(t)\ < Me at 
e at h(t ) 

e at t n , n — 1,2,--- 
e Qt sin (art) 
e at cos (cat) 

f{t — a)h(t — a), a > 0 
with \f(t)\ < Me at 


F{s — a), s > a + a 


_i 

S—O' 


—, 5 > Q: 

— O! ’ 

Ti! ^ 

(s-a) n+1 ’ S ^ a 

t - w i o , s > a 

(s—a) z +ui^ 7 

7— S yf i > ; S > a 
(s—a) z -\-uj z 1 

e~ as F{s), s > a 
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f(t) 

F(s) (continued) 

hit — a), a > 0 

p r\ 

s > S>0 

tf(t) 

-F’is) 

J- sin ut 

2 uj 

{s‘ 2 +uj‘ 2 Y ) s > 0 

^3 [sin ut — ut cos ut] 

( s 2 +aJ 2)2 , S > 0 

f'(t), with f[t) continuous 

sF(s) - fi 0 ) 

and \f(t)\ < Me at 

s > max{a, 0} + 1 

with f\t ) continuous 

S 2 F( S )- S /(0)-/'(0) 

and |/"(t)| < Me at 

s > max{a, 0} + 1 

/(”)(£), with / ( ' n-1 ^(£) continuous 

s n Fis ) — s n ~ l fi 0) — • • • 

and /^(t) < Me at 

-sf n ~ 2 \0) -/ (n - 1} (0) 
s > max{a, 0} + 1 

f* fiu)du, with \f(t)\ < Me at 

s > max{a, 0} + 1 


Table C 
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Practice Problems 

Problem 44.1 

Use Table C to find C[ 2e t + 5]. 

Problem 44.2 

Use Table C to find C[e 3t ~ 3 h(t — 1)]. 

Problem 44.3 

Use Table C to find £ [sin 2 cut]. 

Problem 44.4 

Use Table C to find £[sin3tcos3f]. 

Problem 44.5 

Use Table C to find £[e 2t cos3f]. 

Problem 44.6 

Use Table C to find C[e At (t 2 + 3f + 5)]. 

Problem 44.7 

Use Table C to find £~ 1 [ s3 1 < 0 95 + j^g]. 

Problem 44.8 

Use Table C to find £~ 1 [ ^_ 5 3 ^ 4 ], 

Problem 44.9 

Use Table C to find £~ 1 [|^]. 

Problem 44.10 

Use Table C to find £ -1 [ e 

Problem 44.11 

Graph the function f(t) = h(t — 1) + h(t — 3) for t > 0, where h(t) is the 
Heaviside step function, and use Table C to find C[f{t)]. 

Problem 44.12 

Graph the function /(f) = t[h(t — 1) — h(t — 3)] for t > 0, where h(t) is the 
Heaviside step function, and use Table C to find C[f{t)]. 
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Problem 44.13 

Graph the function /(f) = 3 [h(t — 1) — h(t — 4)] for f > 0, where h(t) is the 
Heaviside step function, and use Table C to find C[f(t)]. 

Problem 44.14 

Graph the function /(f) = |2 — t\[h(t — 1) — h(t — 3)] for f > 0, where h(t) is 
the Heaviside step function, and use Table C to find C[f(t)]. 

Problem 44.15 

Graph the function /(f) = h(2 — t) for f > 0, where h(t) is the Heaviside step 
function, and use Table C to find C[f{t)]. 

Problem 44.16 

Graph the function /(f) = h(t — 1) + h( 4 — f) for f > 0, where h(t) is the 
Heaviside step function, and use Table C to find £[/(f)]. 

Problem 44.17 

The graph of /(f) is given below. Represent /(f) as a combination of Heav¬ 
iside step functions, and use Table C to calculate the Laplace transform of 

f(t). 



Problem 44.18 

The graph of /(f) is given below. Represent /(f) as a combination of Heav¬ 
iside step functions, and use Table C to calculate the Laplace transform of 
/(*)■ 


0 1 


2 



*> 
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Problem 44.19 

Using the partial fraction decomposition find £ _1 


Problem 44.20 

Using the partial fraction decomposition find L~ l 


24e~ 5a 
s 2 —9 


Problem 44.21 

Use Laplace transform technique to solve the initial value problem 


where 


y' + 4y = g(t), ?/(0) = 2 


f 0, 0 < t < 1 
g{t) = \ 12, 1 < t < 3 
I 0, t > 3 


Problem 44.22 

Use Laplace transform technique to solve the initial value problem 
y"-Ay = e 3t , 2/(0) = 0, 2/'(0) = 0. 


Problem 44.23 

Obtain the Laplace transform of the function f 2 tf(\)d\ in terms of £[/(£)] 
F(s) given that J Q 2 f(X)d\ = 3. 
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45 The Laplace Transform and the Method 
of Partial Fractions 


In the last example of the previous section we encountered the equation 

s 2 + 1 

s 2 (s 2 — 4s + 9) 

We would like to find an explicit expression for y(t). This can be done using 
the method of partial fractions which is the topic of this section. According 
to this method, finding ZD 1 , where N(s) and D(s) are polynomials, 

require decomposing the rational function into a sum of simpler expressions 
whose inverse Laplace transform can be recognized from a table of Laplace 
transform pairs. 

The method of integration by partial fractions is a technique for integrating 
rational functions, i.e. functions of the form 


v(t) = C ' 


R(s) 


N(s) 

D(s) 


where N(s) and D(s) are polynomials. 

The idea consists of writing the rational function as a sum of simpler frac¬ 
tions called partial fractions. This can be done in the following way: 


Step 1. Use long division to find two polynomials r(s) and q(s) such that 

N{s) = , w + r(s) 


D(s) ' D( S y 

Note that if the degree of N(s) is smaller than that of D(s) then q(s) = 0 
and r(s) = N(s). 


Step 2. Write D(s) as a product of factors of the form (as + b) n or (as 2 + 
bs + c) n where as 2 + bs + c is irreducible, i.e. as 2 + bs + c = 0 has no real zeros. 


Step 3. Decompose into a sum of partial fractions in the following 
way: 

(1) For each factor of the form (s — a) k write 


Ai A 2 

s — a (s — a) 2 


A-k 


+ ••• + 


(s-a) fc ’ 



where the numbers Ai, A 2 , ■ ■ ■ , A k are to be determined. 

(2) For each factor of the form (as 2 + bs + c) k write 

B\S + Ci B 2 S + C*2 5fcS + Cfc 

as 2 + 6s + c (as 2 + bs + c) 2 (as 2 + bs + c) fc ’ 

where the numbers B\, B 2 , ■ ■ • , Bp. and C\, C 2 , ■ ■ ■ ,Ck are to be determined. 

Step 4- Multiply both sides by D(s) and simplify. This leads to an ex¬ 
pression of the form 

r(s) = a polynomial whose coefficients are combinations of AyBj, and Q. 

Finally, we find the constants, Ai, Bi, and by equating the coefficients of 
like powers of s on both sides of the last equation. 

Example 45.1 

Decompose into partial fractions R(s) = 

Solution. 

Step 1. ^+^±2 = s + 1 + A=t- 
Step 2. s 2 *-l = (s-l)(s + l). 

Step 3. (8+ ;^_i) — i+i + i^T- 

Step 4. Multiply both sides of the last equation by (s — l)(s + 1) to obtain 

s + 3 = A(s — 1) + B(s + 1). 

Expand the right hand side, collect terms with the same power of s, and 
identify coefficients of the polynomials obtained on both sides: 

s + 3 = (A + B)s + (B - A). 

Hence, A + B = 1 and B — A = 3. Adding these two equations gives B — 2. 
Thus, A = — 1 and so 

s 3 + s 2 + 2 , 1 1 , 2 

- 2 - 1 — — s + 1 -- -I- 

s z — 1 s + 1 s — 1 

Now, after decomposing the rational function into a sum of partial fractions 
all we need to do is to ford the Laplace transform of expressions of the form 

A „„ Bs-\-C 
(s—a) n ( as 2 -\-bs-\-c) n * 


29 



Example 45.2 

Find C- 1 


Solution. 

We write 

1 A B 


s{s- 3) s s — 3 

Multiply both sides by s(s — 3) and simplify to obtain 

1 = H(s — 3) + Bs 


or 

1 = (A + B)s-3A. 


Now equating the coefficients of like powers of s to obtain —3 A = 1 and 
A + B = 0. Solving for A and B we find A = — | and B — Thus, 


c - 1 


i 

s(s— 3) 


- W [;] + W 1 [ik] 

— | h(t) + ^e 3 *, t > 0 


where h(t) is the Heaviside unit step function ■ 


Example 45.3 

FM C-‘ [|f±J] . 

Solution. 

We factor the denominator and split the integrand into partial fractions: 


3s + 6 _ A B 
s(s + 3) s s + 3 

Multiplying both sides by s(s + 3) to obtain 


3s + 6 = H(s + 3) + Bs 
= (A + B)s + 3H 


Equating the coefficients of like powers of x to obtain 3 A = 6 and A + B = 3. 
Thus, A — 2 and B = 1. Finally, 


3s + 6 
s 2 + 3s 



+ /T 1 


1 

s + 3 


2 h(t) + e 3t , t> 0. ■ 
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Example 45.4 


Find C~ 


s(s+l ) 2 


Solution. 

We factor the denominator and split the rational function into partial frac¬ 
tions: 

s 2 + 1 _ A B C 

s(s + l) 2 s + s + 1 + (s + 1) 2 ‘ 

Multiplying both sides by s(s + l) 2 and simplifying to obtain 


s 2 + 1 = kL(s + l) 2 + Bs(s + 1) + Cs 

= (A + B)s 2 + (2 A + B + C)s + A. 


Equating coefficients of like powers of s we find A = 1,2 A + B + C = 0 
and A + B = 1. Thus, B — 0 and C = —2. Now hireling the inverse Laplace 
transform to obtain 


C - 1 


s 2 + l 
s(s + l) 2 



-2 Cr 1 


1 

(s + l) 2 _ 


h(t) — 2 te l , t > 0. ■ 


Example 45.5 

Use Laplace transform to solve the initial value problem 


y" + 3y' + 2y = e- t , y(0) = j/(0) = 0. 


Solution. 

By the linearity property of the Laplace transform we can write 

C[y"} + 3 C[y'] + 2 C[y\ = C^]. 

Now since 

C[y"] = s 2 C[y} - sy{ 0) - y'{ 0) = s 2 E(s) 

C[y'} = sY(s) - y( 0) = sE (s) 

= * 
where C[y) = E(s), we obtain 

s 2 E(s) + 3sE(s) + 2Y(s) = 

5 + 1 
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Rearranging gives 


1 


Thus, 


and 


(s 2 + 3s + 2)y (s) — 


Y(s) = 


s + 1 


y{t) = £ 1 


(s + l)(s 2 + 3s + 2)' 
1 


L(s + 1)(s 2 + 3s + 2) 
Using the method of partial fractions we can write 


Thus, 

y(t) = £~ l 


1 1 1,1 

(s + l)(s 2 + 3s + 2) ““ s + 2 ~ s T 1 + (s + 1) 


1 

s + 2 


-/r 1 


i 

s + 1 


+C - 1 


1 

(s + 1) 2 


= e 


- 2 1 




+te \ t> Oi 
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Practice Problems 


In Problems 45.1 - 45.4, give the form of the partial fraction expansion for 
F(s). You need not evaluate the constants in the expansion. However, if the 
denominator has an irreducible quadratic expression then use the completing 
the square process to write it as the sum/difference of two squares. 


Problem 45.1 


Problem 45.2 


Problem 45.3 


Problem 45.4 


Problem_45.5 

Find C~ l i 


Problem 45.6 

Find £-‘ USd 


Problem 45.7 

Find C- 1 [^±1 


Problem 45.8 


Find 


s 2 +6s+8 

s 4 +8s 2 +16 


F(s) = 


s 3 + 3s + 1 
(s — l) 3 (s — 2) 2 ' 


F(s) = 


s 2 + 5s - 3 
(s 2 + 16) (s — 2) 


F(s) = 


s 3 - 1 


(s 2 + l) 2 (s + 4) 2 


F(s) = 


s 4 + 5s 2 + 2s - 9 
(s 2 + 8s + 17)(s-2) 5 
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Problem 45.9 

Use Laplace transform to solve the initial value problem 

y' + 2y — 26 sin 3 t, y( 0) = 3. 

Problem 45.10 

Use Laplace transform to solve the initial value problem 

y' + 2y = 4 1, y( 0) = 3. 

Problem 45.11 

Use Laplace transform to solve the initial value problem 

y" + 3 y' + 2 y = 6e~\ y( 0) = 1, V(0) = 2. 

Problem 45.12 

Use Laplace transform to solve the initial value problem 
y" + 4 y = cos 2 1, y{ 0) = 1, y'{ 0) = 1. 

Problem 45.13 

Lise Laplace transform to solve the initial value problem 
y"-2y'+ y = e 2t , y( 0) = 0, y'(0) = 0. 


Problem 45.14 

Lise Laplace transform to solve the initial value problem 
y" + 9y = g(t), j/(0) = 1, j/(0) = 0 


where 



0 < t < n 

71 < t < OO 


Problem 45.15 

Determine the constants a, (3, y 0 , and y' 0 so that Y(s) = is the Laplace 

transform of the solution to the initial value problem 


y" + ay' + f3y = 0, y(0) = y 0 , y\ 0) = y' 0 . 


Problem 45.16 

Determine the constants a,/3,y 0 , and y' Q so that U(s) = (s + 1)2 is the Laplace 
transform of the solution to the initial value problem 

V" + ay 1 + /3y = 0, y{ 0) = y 0 , y\ 0) = y' 0 . 
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46 Laplace Transforms of Periodic Functions 


In many applications, the nonhomogeneous term in a linear differential equa¬ 
tion is a periodic function. In this section, we derive a formula for the Laplace 
transform of such periodic functions. 

Recall that a function f(t ) is said to be T—periodic if we have f(t+T) = f(t) 
whenever t and t + T are in the domain of fit). For example, the sine and 
cosine functions are 2n —periodic whereas the tangent and cotangent func¬ 
tions are n— periodic. 

If f{t) is T —periodic for t > 0 then we define the function 

fit), 0 <t<T 
0, t > T 

The Laplace transform of this function is then 

poo pT 

C[f T (t)} = / f T (t)e~ st dt = / fit)e~ st dt. 

Jo Jo 

The Laplace transform of a T —periodic function is given next. 



Theorem 46.1 

If fit) is a T— periodic, piecewise continuous fucntion for t > 0 then 


£[/(*)] 


Qfriff] 

1 — e ~ sT ’ 


s > 0. 


Proof. 

Since fit) is piecewise continuous, it is bounded on the interval 0 < t < T. 
By periodicity, fit) is bounded for t > 0. Hence, it has an exponential order 
at infinity. By Theorem 43.2, £[/(£)] exists for s > 0. Thus, 

/»oo 00 f*T 

£[/(£)] = / fit)e~ st dt = ' s ^ / fx{t — nT)hit — nT)e~ st dt, 

Jo n=0 Jo 

where the last sum is the result of decomposing the improper integral into a 
sum of integrals over the constituent periods. 

By the Second Shifting Theorem (i.e. Theorem 44.2) we have 

C[f T {t - nT)h{t - nT )] = e~ nTs C[f T {f)}, s > 0 
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Hence, 

oo / oo 

mO] = ^e-” T ‘£[/r(i)] = £[/r(<)] E'"" 7 '* 

n =0 \n=0 

Since s > 0, it follows that 0 < e~ nTs < 1 so that the series e~ nTs is a 
convergent geoemetric series with limit 1 _ ( 1 _ sT . Therefore, 



mt)\ 


Writ)] 

l-e~ sT ' 


s > 0 ■ 


Example 46.1 

Determine the Laplace transform of the function 

f 1, 0 <t < f 

f(t)={ “ f(t + T) = f(t),t> 0. 

[ 0, | < t < T 

Solution. 

The graph of f(t ) is shown in Figure 46.1. 



By Theorem 46.1, 


£[/(*)] 


fo 2 e St(]t 


1 - e- sT ’ 


Evaluating this last integral, we find 


s > 0. 


mt)] 


sT 

l-e~^r 


o-sT 


1 

—- W7> s > 0 ■ 

s(l + e 2 ) 


Example 46.2 

Find the Laplace transform of the sawtooth curve shown in Figure 46.2 
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M) 



The given function is periodic of period b. For the first period the function 
is defined by 

fb(t) = jt[h(t) — h(t — b)]. 


So we have 


But 


Hence, 


Finally, 


c[h(t)} = £[ff(fc(f) - h(t - 6))] 

-f ±C\h{t) - hit - b)\ 

C[h(t) — h(t — b)\ = C[h(t)\ — C[h(t — b)\ 

1 e~ ba ^ 


cm)\ 


a ( 1 bse bs + e bs 
b ^s 2 s 2 




cm)] 

1 - e~ bs 


a 

b 


1 — e bs — bse bs 
s 2 (l — e ~ bs ) 


Example 46.3 

Find C- 1 [4 - 


e~ a 

s(l-e- s ) 


Solution. 

Note first that 

1 e~ s 1 — e~ s — se~ s 

s 2 s(l - e~ s ) s 2 (l - e~ s ) 
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According to the previous example with a = 1 and b = 1 we find that 
C ^2 — s (iLe-s) is the sawtooth function shown in Figure 46.2 ■ 

Linear Time Invariant Systems and the Transfer Function 

The Laplace transform is a powerful technique for analyzing linear time- 
invariant systems such as electrical circuits, harmonic oscillators, optical de¬ 
vices, and mechanical systems, to name just a few. A mathematical model 
described by a linear differential equation with constant coefficients of the 
form 

a n y^ n) + a n -iy {n - 1] H-f a x y' + a 0 y = b m u^ m) + H - f biu + b 0 u 

is called a linear time invariant system. The function y(t) denotes the 
system output and the function u(t) denotes the system input. The system is 
called time-invariant because the parameters of the system are not changing 
over time and an input now will give the same result as the same input later. 
Applying the Laplace transform on the linear differential equation with null 
initial conditions we obtain 

a n s n Y (s)+a n _is n 1 T(s)+- • --|-ao^( s ) — b m s m U[s)-\-b m — \s m 1 f/(s)+- • --\-boU (s). 
The function 

$( s ) = y ( g ) = bmSm + H-f b x s + b 0 

U (5) a n s n T ci n _\S n ^ T • • • T ct'iS T clq 

is called the system transfer function. That is, the transfer function of 
a linear time-invariant system is the ratio of the Laplace transform of its 
output to the Laplace transform of its input. 

Example 46.4 

Consider the mathematical model described by the initial value problem 

my" + 7 y' + ky = f(t), y( 0) = 0, 2 /(°) = 0. 

The coefficients m, 7, and k describe the properties of some physical system, 
and /(f) is the input to the system. The solution y is the output at time f. 
Find the system transfer function. 
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Solution. 

By taking the Laplace transform and using the initial conditions we obtain 

(ms 2 + 7 s + k)Y(s) = F(s). 


Thus, 


*(s) 


Y(s) 
F(s) 


ms 2 + 7 s + k 


( 2 ) 


Parameter Identification 

One of the most useful applications of system transfer functions is for system 
or parameter identification. 


Example 46.5 

Consider a spring-mass system governed by 

my" + 7 V +ky = f(t ), y(0) = 0, y\ 0) = 0. (3) 

Suppose we apply a unit step force fit) = h(t ) to the mass, initially at 
equilibrium, and you observe the system respond as 

. , 1 , 1 _ t . 1 
y(t) = -~ e cost — -e sin t +-. 

What are the physical parameters m, 7 , and k? 


Solution. 

Start with the model (3)) with /(f) = h(t) and take the Laplace transform of 
both sides, then solve to find Y(s) = 777777 ^ 7 /) • Since /(f) = h(t), F(s) = C 
ffence 


4 ( S ) 


Y(s) 

F(s) 


1 

ms 2 + 7 s + k 


On the other hand, for the input /(f) = h(t) the corresponding observed 
output is 

1 1 . 1 

y{t) = ~ g 6 cosf--e sinf + -. 

ffence, 

E(s) = C[—ke~ f cost — |e _t sin f + |] 

1 s+l _ 1 _ 1 , j_" 

2 (s+l) 2 +l 2 (s+l) 2 +l " 1_ 2s 

s(s 2 +2s+2) 
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Thus, 


4>w 


IM 

F(s) 


By comparison we conclude that m — 1 


1 

s 2 + 2s + 2' 

,7 = 2, and k — 2 ■ 
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Practice Problems 


Problem 46.1 

Find the Laplace transform of the periodic function whose graph is shown. 



Problem 46.2 

Find the Laplace transform of the periodic function whose graph is shown. 



Problem 46.3 

Find the Laplace transform of the periodic function whose graph is shown. 
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Problem 46.4 

Find the Laplace transform of the periodic function whose graph is shown. 


2 



2 3 4 5 6 ^8 


Problem 46.5 

State the period of the function /(f) and find its Laplace transform where 


sin f, 0 < t < 7T 


m 


f(t + 2tt) = fit), t > 0. 


0, 7T < t < 2 tt 


Problem 46.6 

State the period of the function f(t) = 1 — e _t , 0 < t < 2, f(t + 2) = /(f), 
and find its Laplace transform. 

Problem 46.7 

Using Example 44.3 find 


Problem 46.8 

An object having mass m is initially at rest on a frictionless horizontal surface. 
At time f = 0, a periodic force is applied horizontally to the object, causing 
it to move in the positive x-direction. The force, in newtons, is given by 



The initial value problem for the horizontal position, x{t), of the object is 


mx”it ) = /(f), x(0) = x'(0) = 0. 
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(a) Use Laplace transforms to determine the velocity, v(t) = x'(t), and the 
position, x(t), of the object. 

(b) Let m — 1 kg, /o — 1 N, and T — 1 sec. What is the velocity, v, and 
position, x, of the object at t — 1.25 sec? 

Problem 46.9 

Consider the initial value problem 

ay" + by +cy = f(t), y(0) = y'(0) = 0, f>0 

Suppose that the transfer function of this system is given by <3>(s) = 9 A . 2 + 1 5s+ 9 • 

(a) What are the constants a, b, and c? 

(b) If f{t) = e~ l , determine F(s), Y(s), and y(t). 

Problem 46.10 

Consider the initial value problem 

ay" + by’ + cy = f(t), y(0) = y'(0) = 0, t> 0 

Suppose that an input /(f) = f, when applied to the above system produces 
the output y(t) = 2(e _t — 1) + t(e + 1), t > 0. 

(a) What is the system transfer function? 

(b) What will be the output if the Heaviside unit step function /(f) = h(t) 
is applied to the system? 


Problem 46.11 

Consider the initial value problem 

y" + y' + y = f(t), y(0) = y'(0) = 0, 


where 


[ 1, 0 < f < 1 

f(t) = { f(t + 2) 

1-1, 1 < f < 2 


(a) Determine the system transfer function <3>(s). 

(b) Determine Y(s). 


m 


Problem 46.12 

Consider the initial value problem 

y'" - % = e t + f, 1 /( 0 ) = y\ 0) = y"i 0) = 0. 

(a) Determine the system transfer function $( 5 ). 

(b) Determine Y(s). 
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Problem 46.13 

Consider the initial value problem 

y" + by' + cy= h(t ), y{0 ) = y 0 , y\0 ) = y’ 0 , t > 0. 

Suppose that C[y(t)] = Y(s) = . Determine the constants b, c, yo, 

and y' 0 . 
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47 Convolution Integrals 

We start this section with the following problem. 


Example 47.1 

A spring-mass system with a forcing function f(t) is modeled by the following 
initial-value problem 

mx" + kx = f(t), x(0) = xq, x'(0 ) = x' 0 . 

Find solution to this initial value problem using the Laplace transform method. 


Solution. 

Apply Laplace transform to both sides of the equation to obtain 
ms 2 X(s) — msxo — mx' 0 + kX(s ) = F(s). 


Solving the above algebraic equation for X (s) we find 


X(s) 


F(s) , msxo , mx' 0 

ms 2 +k * ms 2 +k ~ ms 2 +k 


1 F i s ) i sxp , 
m S 2 + A ^ s 2 + F ^ s 2 + A - 

m m m 


Apply the inverse Laplace transform to obtain 


x(t) = 




-i 


s 2 + F 


£-'[ *(»)] 

+ Xq C 1 { _ 2 f fe | + x' 0 £ 


-1 


- 




= 1 { F ( s )-^x}+ a; oCOs( v /J)t + x / OA /fsin( v /A)t 

Finding £ _1 |f(s) • s2 |± |,i e., the inverse Laplace transform of a product, 
requires the use of the concept of convolution, a topic we discuss in this 
section ■ 

Convolution integrals are useful when hireling the inverse Laplace transform 
of products H(s) = F(s)G(s). They are defined as follows: The convolution 
of two scalar piecewise continuous functions /(f) and g(t ) defined for t > 0 
is the integral 

( f*g)(t) = [ f(t- s)g(s)ds. 

Jo 
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Example 47.2 

Find f * g where f(t) = e _t and g(t ) = sin t. 

Solution. 

Using integration by parts twice we arrive at 

(/ * 9 ) (t) — fo e~^~ s ^ sin sds 

= \ [ e - d-s) (sin s — cos s)] * 

= ^ + \ (sin t — cos t) ■ 

Graphical Interpretation of Convolution Operation 

For the convolution 

(. f*9)(t) = [ f(t ~ s)g(s)ds 
Jo 

we perform the following: 

Step 1. Given the graphs of /(s) and g(s). (Figure 47.1(a) and (b)) 

Step 2. Time reverse /(—s). (See Figure 47.1(c)) 

Step 3. Shift /(— s) right by an amount t to get f(t — s). (See Figure 47.1(d)) 
Step 4. Determine the product f(t — s)g(s). (See Figure 47.1(e)) 

Step 5. Determine the area under the graph of f(t — s)g(s) between 0 and t. 
(See Figure 47.1(e)) 




f(s) f(-s) 


Vg(s) 


0 s 0 

s t 

s 


(a) 

(b) 

(c) 


f(t-s) ‘ 

|/ vw 


0 t s d 

t 



«») 

(e) 



Figure 47.1 


Next, we state several properties of convolution product, which resemble 
those of ordinary product. 
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Theorem 47.1 

Let f(t),g(t), and k(t) be three piecewise continuous scalar functions defined 
for t > 0 and C\ and C 2 are arbitrary constants. Then 

(i) f * g = g * f (Commutative Law) 

(ii) (f * g) * k = f * (g * k) (Associative Law) 

(iii) / * (cig + c 2 k ) = Cif * g + c 2 f * k (Distributive Law) 

Proof. 

(i) Using the change of variables r = t — s we find 

(. f*g)(t ) = Jof(t-s)g{s)ds 

~ It /( r )^( t - T ) dr 

= Io9(t-r)f(r)dr = (g * f)(t) 

(ii) By definition, we have 

[(/* g) *k)}(t) = f*(f*g)(t-u)k(u)du 


= fo fo U f(t - u - w)g(w)k(u)dw 


du 


For the integral in the bracket, make change of variable w = s — u. We have 

ft r rt 


[( f*g ) * k)](t) = 


fit — s)g(s — u)k{u)ds 


du. 


Jo L Ju 

This multiple integral is carried over the region 

{(s, u) : 0 < u < s < t} 

as depicted by shaded region in the following graph. 

u 



s 

Figure 47.2 
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Changing the order of integration, we have 


[(/ *g) * k)](t) = /o [/o' f(t- s)g(s - u)k(u)du] ds 

= Jof( t ~ s )( 9 *k){s)ds 

[f*(g*k)](t ) 

(iii) We have 

(f * (cig + c 2 k))(t) = Jo f(t ~ s)(cig(s) + c 2 k(s))ds 


= c i fo f(t ~ s)g(s)ds + c 2 f* f(t - s)k(s)ds 


= Ci(f * g)(t) + c 2 (f * k)(t) m 

Example 47.3 

Express the solution to the initial value problem y' + ay = g(t), y( 0) = y 0 
in terms of a convolution integral. 


Solution. 

Solving this initial value problem by the method of integrating factor we find 

y(t) = e~ at y 0 + [ e~ a{t ~ s) g(s)ds = e~ Qt y 0 + e~ at * g(t) u 

Jo 

Example 47.4 

If f(f) is an m x n matrix function and g(t) is an n X p matrix function then 
we define 

(f*g )(t)=[ f(t - s)g(s)ds, t > 0. 

Jo 

Express the solution to the initial value problem y' = Ay + g (t), y(0) = y 0 
in terms of a convolution integral. 


Solution. 

The unique solution is given by 

y(t) = e ,A y„ + /V ( ‘- S) g (s)ds = e‘ A y„ + e ,A * g(t) 1 
Jo 

The following theorem, known as the Convolution Theorem, provides a way 
for finding the Laplace transform of a convolution integral and also finding 
the inverse Laplace transform of a product. 
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Theorem 47.2 

If f{t) and g(t ) are piecewise continuous for t > 0, and of exponential order 
at infinity then 

£[(/ * £)(*)] = £[/(f)]£[tf(f)] = F ( S ) G ( S )- 

Thus, (/^)(t) = £- 1 [F( S )G( S )]. 


Proof. 

First we show that f * g has a Laplace transform. From the hypotheses we 
have that \f(t)\ < Mie ait for t > C\ and | g(t)\ < M 2 e a2t for t > C 2 . Let 
M = M\M 2 and C = C\ + C 2 . Then for t > C we have 


\(f*9)(t)\ 


Iof(t-s)g{s)ds < \f(t 


s)||^(s)|ds 


< M\M 2 fg e aid s )e a2S ds 


Mte ait , 
M ea2t ~ eait 

a. 2 —ai ’ 


(I | — cl, 2 
Q'l ~f~ 0>2 


This shows that / * g is of exponential order at infinity. Since / and g are 
piecewise continuous, the first fundamental theorem of calculus implies that 
f * g is also piecewise continuous. Hence, f * g has a Laplace transform. 
Next, we have 


£[(f*9)(t)] = / 0 °° e (fo fit ~ T )g( T ) dT ) dt 


= /t = 0 f*=o e st f(t-r)g(r)dTdt 

Note that the region of integration is an infinite triangular region and the 
integration is done vertically in that region. Integration horizontally we find 

/»oo POO 

£[(f*9)(t)] = / / e~ st f(t — T)g(r)dtdT. 

J T=0 J t—T 

We next introduce the change of variables f3 = t—T. The region of integration 
becomes r > 0, i > 0. In this case, we have 

£[(f*9)(t)} = IZoI^o e ~ s(/3+T) f^d(T)dTd/3 

= (fZo e ~ ST g{ T )dr) (fZo e ~ Sf} f(P) d P) 

G(s)F(s) = F(s)G(s ) ■ 
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Example 47.5 

Use the convolution theorem to find the inverse Laplace transform of 


H(s) = 


(s 2 + a 2 ) 2 ' 


Solution. 

Note that 


H(s) = 


s 2 + a 2 J \ s 2 + a 2 


So, in this case we have, F(s) = G(s ) = s2 + a 2 so that /(f) = g(t ) = 1 sin (at). 
Thus, 


1 [ f 1 

(/ * g)(t) = —t / sin (at — as) sin (as)ds = —Usin (at) — at cos (at)) ■ 

a- J 0 2 a 6 

Convolution integrals are useful in solving initial value problems with forcing 
functions. 


Example 47.6 

Solve the initial value problem 

4y" + y = g(t), y(0) = 3, y'(0) = -7 

Solution. 

Take the Laplace transform of all the terms and plug in the initial conditions 
to obtain 

4(s 2 Y'(s) - 3s + 7) + Y'(s) = G(s) 


or 


(4s 2 + l)y(s) - 12s + 28 = G(s). 
Solving for Y(s) we find 

Y(8) = 


12s—28 , G(s) 

4U 2 + U + 4( S 2+U 


ffence, 


_ 3s _ 7 ( 2 ) 1 ( 2 ) 

* 2 +((§) 2 + 4 T ’ s 2 + G ) 2 

y(t) =3 cos (^j -7sin^0 + \J Q (£) g(t - s)ds. 


So, once we decide on a g(t) all we need to do is to evaluate the integral and 
we’ll have the solution ■ 
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Practice Problems 


Problem 47.1 

Consider the functions f(t) = g(t ) = h(t), t > 0 where h(t) is the Heaviside 
unit step function. Compute f * g in two different ways. 

(a) By directly evaluating the integral. 

(b) By computing £~ 1 [F(s)G(s)] where F(s ) = £[/(£)] and G(s ) = C[g(t)]. 

Problem 47.2 

Consider the functions f(t) = e l and g{t ) = e~ 2t , t > 0. Compute f * g in 
two different ways. 

(a) By directly evaluating the integral. 

(b) By computing C~ 1 [F(s)G(s)] where F(s ) = £[/(£)] and G(s) = C[g(t)]. 

Problem 47.3 

Consider the functions f(t) = sin t and g(t) = cos t, t > 0. Compute f * g in 
two different ways. 

(a) By directly evaluating the integral. 

(b) By computing £~ 1 [F(s)G(s)] where F(s) = C[f(t)] and G(s) = C[g(t)]. 


Problem 47.4 

Use Laplace transform to comput the convolution P * y, where | bfP(t) = 


h(t) 

0 


e 

t 


and y (t) — 


h(t ) 


Problem 47.5 

Compute and graph f * g where f(t) = h(t) and g(t) = t[h(t ) — h(t — 2)]. 

Problem 47.6 

Compute and graph f * g where f(t) = h(t) — h(t — 1) and g(t) = h(t — 1) — 
2h(t — 2)]. 


Problem 47.7 

Compute t*t*t. 


Problem 47.8 

Compute h(t) * e~ l * e~ 2t . 

Problem 47.9 

Compute t * e _t * e l . 
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Problem 47.10 


n functions 


/ ‘ ^N 

Suppose it is known that hit) * h(t ) * • • • * h(t ) = Ct 8 . Determine the con¬ 
stants C and the poisitive integer n. 

Problem 47.11 

Use Laplace transform to solve for y(t) : 

sin (t — X)y(X)d\ = t 2 . 



Problem 47.12 

Use Laplace transform to solve for y{t) : 

y{t) - [ e (t ~ X) y(X)dX = t. 
Jo 


Problem 47.13 

Use Laplace transform to solve for y(t) : 

t*y(t ) = t 2 (l - e -< ). 


Problem 47.14 

Use Laplace transform to solve for y(t) : 


y = h{t) * y, y(0) 


1 

2 


Problem 47.15 

Solve the following initial value problem. 

y'-y= I (t-X)e x dX, y( 0) =-1. 
Jo 
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48 The Dirac Delta Function and Impulse Re¬ 
sponse 

In applications, we are often encountered with linear systems, originally at 
rest, excited by a sudden large force (such as a large applied voltage to an 
electrical network) over a very short time frame. In this case, the output 
corresponding to this sudden force is referred to as the ”impulse response”. 
Mathematically, an impulse can be modeled by an initial value problem with 
a special type of function known as the Dirac delta function as the external 
force, i.e., the nonhomogeneous term. To solve such IVP requires finding the 
Laplace transform of the delta function which is the main topic of this section. 

An Example of Impulse Response 

Consider a spring-mass system with a time-dependent force /(f) applied to 
the mass. The situation is modeled by the second-order differential equation 

my" + ~/y' + ky = f(t) (4) 

where f is time and y(t ) is the displacement of the mass from equilibrium. 
Now suppose that for f < 0 the mass is at rest in its equilibrium position, so 
y(0) = y'(0) = 0. Hence, the situation is modeled by the initial value problem 

my"+ 7 P + ky = f(t), y(0) = 0, y'(0) = 0. (5) 

Solving this equation by the method of variation of parameters one finds the 
unique solution 



4>(t — s)f(s)ds 


( 6 ) 


where 


</>(f) 


e ( 7/2m)t g j n 




m 



7 


2 


4m 2 


Next, we consider the problem of strucking the mass by an ”instantaneous” 
hammer blow at f = 0. This situation actually occurs frequently in practice-a 
system sustains a forceful, almost-instantaneous input. Our goal is to model 
the situation mathematically and determine how the system will respond. 
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In the above situation we might describe f(t) as a large constant force applied 
on a very small time interval. Such a model leads to the forcing function 

at) = { -e’ 

ew | 0, otherwise 

where e is a small positive real number. When e is close to zero the applied 
force is very large during the time interval 0 < t < e and zero afterwards. A 
possible graph of f e (t) is given in Figure 48.1 


\le 


£ 

Figure 48.1 


In this case it’s easy to see that for any choice of e we have 



1 


and 


lim f e (t) = 0, t ± 0, lim f e ( 0) = oo. (7) 

e—>-0+ £^-0+ 

Our ultimate interest is the behavior of the solution to equation (4) with 
forcing function f e (t) in the limit e —> 0 + . That is, what happens to the 
system output as we make the applied force progressively ’’sharper” and 
’’stronger?”. 

Let y e (t) be the solution to equation (4) with /(f) = f e (t). Then the unique 
solution is given by 

Ve{t) = [ s)f £ (s)ds. 

■Jo 
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For t > e the last equation becomes 

Ve(t) = ~ I s)ds. 

e Jo 

Since 4>(t) is continuous for all t > 0 we can apply the mean value theorem 
for integrals and write 

Ve(t) = 

for some 0 < if < e. Letting e —> 0 + and using the continuity of <f> we find 

y(t) = lim y e (t) = f(t). 

€— > 0 + 

We call y(t) the impulse response of the linear system. 

The Dirac Delta Function 

The problem with the integral 


</>(t - s)f e (s)ds 

is that lim e _ 5 . 0 + / e (0) is undefined. So it makes sense to ask the question of 
whether we can find a function 5(t) such that 

lim € ^ 0 + y e (t) = lim £ _ >0 + 4>{t — s)f e {s)ds 

= (j)(t — s)5(s)ds 

where the role of 5(t) would be to evaluate the integrand at s — 0. Note that 
because of Fig 48.1 and (7), we cannot interchange the opeartions of limit 
and integration in the above limit process. Such a function 6 exist in the 
theory of distributions and can be defined as follows: 

If f{t) is continuous in a < t < b then we define the function 5(t) by the 
integral equation 




f(t)6(t - t Q )dt 


lim f f(t)f e (t - t 0 )dt. 


The object 5{t) on the left is called the Dirac Delta function, or just the 
delta function for short. 
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Finding the Impulse Function Using Laplace Transform 

For e > 0 we can solve the initial value problem (5) using Laplace transforms. 
To do this we need to compute the Laplace transform of f e (t), given by the 
integral 


r°° 1 r e 1 _ P -es 

C[f £ (t )] = / Mt)e~ st dt = - / e~ st dt = -—. 

Jo e Jo 

Note that by using L’Hopital’s rule we can write 

lim C[fe(t)] = lim ---- = 1 , s > 0 . 

e^o+ v ' e^o+ es 

Now, to find y e (t), we apply the Laplace transform to both sides of equation 
(4) and using the initial conditions we obtain 

1 — e~ es 

ms 2 YAs ) + 'ysYJs) + kYAs ) =--—. 

es 

Solving for F e (s) we find 


Y e (s) 


1 1 - e- £S 

ms 2 + ys + k es 


Letting e —» 0 + we find 


Y{s) 


1 

ms 2 + ys + k 


which is the transfer function of the system. Now inverse transform Y(s) to 
hnd the solution to the initial value problem. That is, 


y(t) = £ 1 ( —y—-— —r) = 

\ms 2 + 7 s + k J 

Now, impulse inputs are usually modeled in terms of delta functions. Thus, 
knowing the Laplace transform of such functions is important when solving 
differential equations. The next theorem finds the Laplace transform of the 
delta function. 
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Theorem 48.1 

With 5(t) defined as above, if a < t 0 < b 

- t 0 )dt = f(t 0 ). 



Proof. 

We have 


la f(t)S(t- to) = linw,+ fa f(t)fe(t ~ t 0 )dt 

= lim € ^ 0+ \ / t * 0+e f(t)dt 

= lim e ^ 0 + \ f{t Q + /3e)e = f(t 0 ) 

where 0 < fd < 1 and the mean-value theorem for integrals has been used ■ 

Remark 48.1 

Since p e (t — t 0 ) = \ for t 0 < t < t 0 + e and 0 otherwise we see that f a — 

a)dt = /(a) and J ^ f(t)5(t — to)dt = 0 for t 0 > b. 

It follows immediately from the above theorem that 

noo 

jC[S(t — t 0 )]= / e~ st 5{t — t 0 )dt = e~ st °, t 0 > 0. 

Jo 

In particular, if to — 0 we find 

cmi = i. 

The following example illustrates the formal use of the delta function. 

Example 48.1 

A spring-mass system with mass 2, damping 4, and spring constant 10 is 
subject to a hammer blow at time t = 0. The blow imparts a total impulse of 
1 to the system, which was initially at rest. Find the response of the system. 

Solution. 

The situation is modeled by the initial value problem 

2y"+ 4y'+ 10y = S(t), y( 0) = 0, y'( 0) = 0. 
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Taking Laplace transform of both sides we find 

2s 2 Y(s) + 4sY» + 10y(s) = 1. 


Solving for Y(s) we find 


Y(s) 


The impulsive response is 


1 

2s 2 + 4s + 10' 


y{t) = £ 1 


2 (s + iy + 2 2 


1 - 9 / 

-e sm 2 1 m 
4 


Example 48.2 

A 16 lb weight is attached to a spring with a spring constant equal to 2 
lb/ft. Neglect damping. The weight is released from rest at 3 ft below the 
equilibrium position. Att = 2n sec, it is struck with a hammer, providing an 
impulse of 4 lb-sec. Determine the displacement function y(t) of the weight. 


Solution. 

This situation is modeled by the initial value problem 
16 

-y" + 2 y = 4 5(t - 2vr), y{ 0) = 3, y\ 0) = 0. 
Apply Laplace transform to both sides to obtain 

s 2 Y» - 3s + 4Y» = 8e -27r L 


Solving for Y (s) we find 


Y{s) = 


3 s 


0 — 27 rs 


s 2 + 4 s 2 + 4 
Now take the inverse Laplace transform to get 

y(t) = £ _1 [Y'( S )] — 3 cos 2 1 + 8 h(t — 2n)f(t — 2i r) 

where 


m = 


i 


= - sin 2 1. 
2 


Hence, 


y(t) = 3 cos 2 1 + 4 h(t — 2 ti) sin 2 (t — 2n) = 3 cos 2 1 + 4 h(t — 2 ti) sin 2 1 


or more explicitly 



3 cos 2 1, 

3 cos 2 1 + 4 sin 2t, 


t < 2n 
t > 2 -k m 
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Practice Problems 

Problem 48.1 

Evaluate 


( a ) Io( 1 + e t )S(t-2)dt. 

(b) J- 2 O- + — 2 )dt. 



cos 2 t 
te _t 


(d)f.y-‘+t) 


8(t)dt. 

5(t + 2 ) 
5(t-l) 
5(t- 3) 


dt. 


Problem 48.2 

Let f{t) be a function defined and continuous on 0 < t < 00 . Determine 

(f *8)(t) = f f(t-s)S(s)ds. 

Jo 

Problem 48.3 

Determine a value of the constant to such that J 0 ' sin 2 \k( t — to)]8(t—^)dt = 


Problem 48.4 

If f^t n S(t — 2)dt = 8, what is the exponent n? 


Problem 48.5 

Sketch the graph of the function git) which is defined by git) = J* J* 5(u — 
1 )duds, 0 < t < 00 . 


Problem 48.6 

The graph of the function g(t) = f* e at S(t — t 0 )dt, 0 < t < 00 is shown. 
Determine the constants a and t 0 . 
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Problem 48.7 

(a) Use the method of integarting factor to solve the initial value problem 

y' -y = h(t), y{ 0) = 0. 

(b) Use the Laplace transform to solve the initial value problem <i>' — <i> — 
5{t), 0(0) = 0. 

(c) Evaluate the convolution (p*h(t ) and compare the resulting function with 
the solution obtained in part (a). 

Problem 48.8 

Solve the initial value problem 

y' + y = 2 + 6(t-l), y( 0) = 0, 0 < t < 6. 

Graph the solution on the indicated interval. 

Problem 48.9 

Solve the initial value problem 

y" = S(t - 1) - 5(t - 3), y(0) = 0, y\ 0) = 0, 0 < t < 6. 

Graph the solution on the indicated interval. 

Problem 48.10 

Solve the initial value problem 

y" - 2 y' = 6(t - 1), y(0) = 1, y\ 0) = 0, 0 < t < 2. 

Graph the solution on the indicated interval. 

Problem 48.11 

Solve the initial value problem 

y" + 2y , + y = 5(t-2), y( 0) = 0, y\ 0) = 1, 0 < t < 6. 

Graph the solution on the indicated interval. 
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49 Solving Systems of Differential Equations 
Using Laplace Transform 


In this section we extend the definition of Laplace transform to matrix-valued 
functions and apply this extension to solving systems of differential equations. 
Let yi(t), y 2 (t), ■ ■ ■ ,y n (t) be members of VS. Consider the vector-valued 
function 


y (*) 


yi(t) 

V2 (t) 


y n (t) 


The Laplace transform of y (t) is 


£[y(t)} = / 0 °° ’y(t) e stdt 


/ 0 °° 2/i(^)e st dt 
y 2 (.t)e- st dt 

r yn(t)e~ st dt 

£[yi(t)] 

zMt)] 

£[y n (t)} 


In a similar way, we define the Laplace transform of an m x n matrix to 
be the m x n matrix consisting of the Laplace transforms of the component 
functions. If the Laplace transform of each component exists then we say 
y(t) is Laplace transformable. 


Example 49.1 

Find the Laplace transform of the vector-valued function 


y(*) = 
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Solution. 

The Laplace transform is 


£[y (*)] 


s > 1 ■ 


1 

S — 1 


The linearity property of the Laplace transform can be used to establish the 
following result. 


Theorem 49.1 

If A is a constant n x n matrix and B is an n x p matrix-valued function 
then 


£[AB(t)] = A£[B(t)]. 


Proof. 

Let A = (a,ij) and B(t) = (%(f)). Then AB(f) = (Y2=i a ikhp)- Hence, 

n n 

£[AB(f)] = dikbkp)\ ^ Q J ikf~'(pkp)\ A£[B(tj] I 

k =1 k =1 


Theorem 42.3 can be extended to vector-valued functions. 

Theorem 49.2 

(a) Suppose that y (t) is continuous for t > 0 and let the components of the 
derivative vector y' be members of VS. Then 

£[y'W] = sC[y(t)] — y(0). 

(b) Let y'(f) be continuous for t > 0, and let the entries of y"(f) be members 
of VS. Then 

£[y"W] = s 2 £[y(t)] - sy(0) - y(0). 

(c) Let the entries of y (t) be members of VS. Then 

£ {i yW *} = ^' 


62 



Proof. 

(a) We have 


c[y'M 

C[y' 2 it) \ 


[ cm)] \ 

s£[yi(t)} — 2 / 1 ( 0 ) 

sC[y 2 {t)} - 1 / 2 ( 0 ) 


sC[y n (t)} - y n ( 0) 
s£[y(t)}-y(0) 


s£[y'(t)]-y>(0) 

= s(s£[y{t)} -y(0)) - y'(0) 

= s 2 £[y(t)]-sy(0)-y'(0) 

(c) We have 

C[y(t)] = sC y(s)rfs| 

so that 

£ {/ yW *} = ^ ■ 

The above two theorems can be used for solving the following initial value 
problem 

y'(t) = Ay+ g(t), y(0)=y 0 , t>0 (8) 

where A is a constant matrix and the components of g (t) are members of 

VS. 

Using the above theorems we can write 

sY(s)-y 0 = AY(s) + G(s) 


(si - A)Y(s) = y 0 + G(s) 


C[y '(*)] = 


(b) We have 


C\y"(t)] = 
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where £[g(t)] = G(s). If s is not an eigenvalue of A then the matrix si — A 
is invertible and in this case we have 

Y(s) = (sI-A)- 1 [y 0 + G(s)]. (9) 

To compute y(t) = £ _1 [Y(s)] we compute the inverse Laplace transform 
of each component of Y(s). We illustrate the above discussion in the next 
example. 


Example 49.2 

Solve the initial value problem 


/ 

y 


1 2 
2 1 



y(o) 


i 

-2 


Solution. 

We have 


and 


Thus, 


(si - A)’ 1 


1 

(s + l)(s - 3) 


s - 1 2 

2 s - 1 


G(s) = 


Y(s) = (sI-A)- 1 [y 0 + G(s)] 


1 

' s - 1 2 


r i 

(s+l)(s-3) 

2 s - 1 


i 

i 

to 

1 

1- 


s 4 -6s 3 +9s 2 -4s+8 
s 2 (s+1)(s-2)(s-3) 
-2s 4 +8s 3 -8s 2 +6s~4 
s 2 (s+l)(s—2)(s—3) 

Using the method of partial fractions we can write 


Yi{s) 

Y 2 (s) 


Therefore 


4 J_ _ 81 , 7 _1_1_L_1 1 

3 s 2 9 s + 3 s+1 3 s—2 9s-3 


2 J_ | 10 1 _ 7 _L_2 _1_1 1 

3s 2 + 9 s 3 s+1 3 s—2 9 s-3 


= U' 1 [Y 2 (-s)] = —ft + f - - |e 2f - 


yi(t) 
V2 (t) 
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fe 3t , t > 0 



Hence, for t > 0 


y{t) = t\ K] + \ J 1 + e-* [ I ] + e 2t I" _| 1 + e 3i T 1 ■ 

3 J L 9 J L3J L 3 J L 9 . 

System Transfer Matrix and the Laplace Transform of e tA 

The vector equation (8) is a linear time invariant system whose Laplace 
input is given by yo + G(s) and the Laplace output Y(s). According to 
(9) the system tranform matrix is given by (si — A)” 1 . We will show that 
this matrix is the Laplace transform of the exponential matrix function e tA . 
Indeed, e tA is the solution to the initial value problem 

$'(t) = A $(t), <L(0) = I, 

where I is the n x n identity matrix and A is a constant n x n matrix. Taking 
Laplace of both sides yields 

s£[$(t)]-I = A£[$(t)]. 

Solving for £[<&(£)] we find 

£[$(£)] = (si-A)’ 1 = C[e tA ]. 
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Practice Problems 


Problem 49.1 

Find £[y(i)] where 

e~' cos 2 1 
0 

t + e l 



Problem 49.2 

Find £[y(i)] where 

y (t) 


Problem 49.3 

Find £ _1 [Y(s)] where 

Y(a) 


Problem 49.4 

Find £ _1 [Y(s)] where 



' 1 -1 2 ' 


- £[f 3 ] 1 

Y (s) = 

2 0 3 


C{e 2t ] 


1-2 1 


£[sinf] 



i 


s 2 +2s+2 

i 

S 2 +S 


Problem 49.5 

Use the Laplace transform to solve the initial value problem 


y 


/ 




o 

o 


Problem 49.6 

Use the Laplace transform to solve the initial value problem 

y 7 = [ 5 ~ 4 

y 3-2 


y, y(o) = 
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Problem 49.7 

Use the Laplace transform to solve the initial value problem 


/ 

y 


1 4 

-1 1 


y + 


o 

3e* 



Problem 49.8 

Use the Laplace transform to solve the initial value problem 


// 


y 


-3 -2 
4 3 


y, y(0) 



o 

l 


Problem 49.9 

Use the Laplace transform to solve the initial value problem 


// 


y 


l 

l 




o 

o 


Problem 49.10 

Use the Laplace transform to solve the initial value problem 


/ 

y = 

' 1 o o' 
0-11 

y + 

e} 

1 

, y(o) = 

" 0 ' 

0 


0 0 2 


—2 1 


0 


Problem 49.11 

The Laplace transform was applied to the initial value problem y' 
y 0 , where y (t) = ^ VlW 


, A is a 2 x 2 constant matrix, and yo 
The following transform domain solution was obtained 

1 


2/2 (t) 
form 

£[y (*)] = Y ( s ) = 


s 2 - 9s + 18 


s — 2 
4 


-1 
s — 7 


2/i,o 

2 / 2,0 


(a) what are the eigenvalues of A? 

(b) Find A. 


Ay, y(0) 

2/i,o 

_ 2 / 2,0 
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50 Solutions to Problems 


Section 43 


Problem 43.1 

Determine whether the integral J 0 °° jr^dt converges. If the integral con¬ 
verges, give its value. 

Solution. 

We have 

f°° 1 f A 1 

/ -- 7,dt = lim / -- -dt ■ 

Jo 1 + t 2 °° Jo 1 + t 2 

, vr 

= lim arctan A = — 

A— >oo 2 

So the integral is convergent ■ 

Problem 43.2 

Determine whether the integral J 0 °° jJ-^dt converges. If the integral con¬ 
verges, give its value. 

Solution. 

We have 

r°° t , i r A 2t , 

/ --77at =- hm / --at 

Jq 1 + t 2 2 -A-»00 1 + t 2 

= - lim In (1 + A 2 ) = 

2 A— »oo 

Hence, the integral is divergent ■ 

Problem 43.3 

Determine whether the integral J 0 °° e~ t cos (e~ t )dt converges. If the integral 
converges, give its value. 


= - lim Tin (1 + f 2 )l n 4 

2a^oo 1 v yj 0 


= hm [arctan f]^ 
A—>oo 



Solution. 

Using the substitution u = e~ l we find 


'o 


e *cos(e t )dt = lim / —cos udu 

y ' A->oo 1 


1 


= lim [— sinu]f = lim [sin 1 — sin (e )] 

A—>• oo A—>-oo 

= sin 1 

Hence, the integral is convergent ■ 

Problem 43.4 

Using the definition, find £[e 3t ], if it exists. If the Laplace transform exists 
then find the domain of F(s). 

Solution. 

We have 

/ A pA 

e 3t e~ st dt= lim / e t{ - 3 ~ s) dt 


A— 


= lim 

A—>00 


= lim 

A—>■ oo 


s-3 


I o 

gt(3—s)' 

3 — s 

e A( 3 -s) 


A —>-oo 


3 — s 3 — s 


, s > 3 


Problem 43.5 

Using the definition, find C[t — 5], if it exists. If the Laplace transform exists 
then find the domain of F(s). 


Solution. 

Using integration by parts we find 

f A 

C[t — 5] = lim / (t — 5 )e st dt = lim 


A—>-oo 


= lim 

A—>-oo 


A—>-oo 


-(t-5)e~ st 


+ - / e st dt 


J o 


— (H — 5)e -sA + 5 


~e~ st ~ 

A \ 

_ 5 2 . 

oj 


s>0 
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Problem 43.6 

Using the definition, find £[e^ _ b ], if it exists. If the Laplace transform 
exists then find the domain of F(s). 


Solution. 

We have 


e{ t-if e - st dt= / e (t “ 1 ) 2 -^f. 


Since lim t _ 5 . 00 (t — l) 2 — st = lirn^oo f 2 ^1 — = oo, for any fixed s 

we can choose a positive C such that (t — l) 2 — st > 0 for t > C. In this case, 
e (t-i) 2 -st > anf ] this implies that / 0 °° e^ t ~ 1 '^~ st dt > J^° dt. The integral on 
the right is divergent so that the integral on the left is also divergent by the 
comparison theorem of improper integrals. Hence, /(f) = does not 

have a Laplace transform ■ 


Problem 43.7 

Using the definition, find £\(t — 2) 2 ], if it exists. If the Laplace transform 
exists then find the domain of F(s). 


Solution. 

We have 


£[(t — 2) 2 ] = lim / (f — 2) 2 e st dt. 

T->oo L 


Using integration by parts with u’ — e st and v — (t — 2) 2 we find 


(f - 2 ) 2 e~ st dt = 


s Jo s Jo 
: 4_ (Tj2fiy + 2^_ 

s S S In 


Thus, 


(f — 2) 2 e st dt = - + - lim [ (f — 2)e st dt. 

S S T— XX) In 


Using by parts with v! — e st and v — t — 2 we find 


(f — 2)e st dt = 


(t - 2)e 
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Letting T —> oo in the above expression we find 


Hence, 


r T 2 1 

lim / (t — 2 )e~ st dt =-b -77, s > 0. 

r^ooy 0 s s 2 



4 

s 



s > 0 ■ 


Problem 43.8 

Using the definition, find £[/(£)], if it exists. If the Laplace transform exists 
then find the domain of F(s). 


m = 


0, 0 < t < 1 

t- 1, t> 1 


Solution. 

We have 


£[/(*)] = J im 

T-> oo 


(1 — l)e dt 


Using integration by parts with u! = e and v — t — 1 we find 


lim / (£ — l)e st dt = lim 

T—>-oo /i T—>• oo 


(t — l)e 


- W e 


— St 


= o - ) S > 0 


Problem 43.9 

Using the definition, find £[/(£)], if it exists. If the Laplace transform exists 
then find the domain of F(s). 


[ 0, 0 < t < 1 

f(t)= l t- 1, 1 < t < 2 
I 0, t > 2. 


Solution. 

We have 


£[/(*)] 


[ (t— l)e Wit = 

(i - l)e- 1 _„1 

9 ^ 

/1 

s s z 


p -2s 1 

-b -^(e _s - e~ 2s ), s / 0 | 

s s z 
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Problem 43.10 

Let n be a positive integer. Using integration by parts establish the reduction 
formula 

r +n p -st r 

/ t n e~ st dt = -— + - / t n ~ 1 e~ st dt, s> 0. 


Solution. 

Let v! = e~ st and v = t n . Then u = — -— and v' = nt n ~ 1 . Hence, 

S 

/ fn -st r 

t n e~ st dt = -— + - t n ~ 1 e~ st dt, s > 0 ■ 


Problem 43.11 

For s > 0 and n a positive integer evaluate the limits 

(a) linif^o t n e~ st (b) lim t ^ 00 t n e~ st 


Solution. 

(a) lini^o t n e~ st = lim^ 0 = f = 0. 

(b) Using L’Hopital’s rule repeatedly we find 
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lim t n e~ st = • • • = lim -- = 0 ■ 

t—>oo t —>oo S n e st 

Problem 43.12 

(a) Use the previous two problems to derive the reduction formula for the 
Laplace transform of f(t) = t n , 

C[t n } = -C[t n ~\ s > 0. 
s 

(b) Calculate C[t k ], for k = 1,2, 3,4, 5. 

(c) Formulate a conjecture as to the Laplace transform of f(t),t n with n a 
positive integer. 


Solution. 

(a) Using the two previous problems we find 


C\t n 1 = lim [ t n e~ st dt = lim 

T—>oc J q T—> oo 


n^—st 


t n e 


J o 


T r T 

+ - / 
s 


: n ~ 1 e~ st dt 


'o 


= - lim 

S T— >oo 


t n - 1 e-* t dt=”C[t n - 1 ], s> 0 


n 
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(b) We have 


m 

£[t 2 ] 

£[t 3 ] 

m 


_ i 

= ~s C ® ^ 

-;<W 



(c) By induction, one can easily show that for n = 


C[t n ] 



s > 0 ■ 


0 , 1 , 2 ,--- 


From a table of integrals, 

I e au sin 0udu =e c 
[ e au cos / 3udu =e c 


, a sin f3u — (3 cos f3u 

a 2 + /3 2 

, a cos (3u + /3 sin f3u 
a 2 + /3 2 


Problem 43.13 

Use the above integrals to find the Laplace transform of f(t) = 
exists. If the Laplace transform exists, give the domain of F(s). 


Solution. 

We have 

£ [cos cat] = lim 

T—> oo 



—s cos cat + ca sin cat 

n 

s 2 + ca 2 

o J 


S 

S 2 + OJ 2 ’ 


Problem 43.14 

Lise the above integrals to find the Laplace transform of f(t) = 
exists. If the Laplace transform exists, give the domain of F(s). 


Solution. 

We have 

£ [sin cat] = lim 

T—>oo 


—st 


-s sin cat + ca cos cat 
s 2 + ca 2 


ca 

s 2 + ca 2 ’ 


coscat, if it 

s > 0 ■ 

sin cat, if it 

s > 0 ■ 
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Problem 43.15 

Use the above integrals to find the Laplace transform of /(f) = cos ca(t — 2), 
if it exists. If the Laplace transform exists, give the domain of F(s). 


Solution. 

Using a trigonometric identity we can write /(f) = cos ca(t — 2) = coscat cos 2ca+ 
sincat sin 2ca. Thus, using the previous two problems we find 


£[cosca(t — 2)] = 


s cos 2ca + ca sin 2ca 
s 2 + ca 2 


s > 0 


Problem 43.16 

Lise the above integrals to find the Laplace transform of f(t) = e 3t sinf, if it 
exists. If the Laplace transform exists, give the domain of F(s). 


Solution. 

We have 



Problem 43.17 

Lise the linearity property of Laplace transform to find £[5e~ 7t + t + 2e 2t ]. 
Find the domain of F(s). 


Solution. 

We have C[e~ 7t \ = s > —7, C[t] = Jj, s > 0, and C[e 2t ] = s > 2. 
Hence, 

£[ 5 e- 7 ' +1 + 2e 2t ] = 5C[e~ 7t } + C[t] + 2 C[e 2t ] = \ s> 2 ■ 

s + 7 s z s — 2 

Problem 43.18 

Consider the function f(t) = taut. 

(a) Is f(t) continuous on 0 < t < oo, discontinuous but piecewise continuous 
on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(t)| < Me at for 0 < t < oo? 
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Solution. 

(a) Since /(f) = tan f = and this function is discontinuous at t = (2 n + 
1)|. Since this function has vertical asymptotes there it is not piecewise 
continuous. 

(b) The graph of the function does not show that it can be bounded by 
exponential functions. Hence, no such numbers a and M ■ 

Problem 43.19 

Consider the function /(f) = t 2 e~ l . 

(a) Is /(f) continuous on 0 < t < oo, discontinuous but piecewise continuous 
on 0 < t < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 

Solution. 

(a) Since f 2 and e _t are continuous everywhere, /(f) = f 2 e _t is continuous on 
0 < f < oo. 

(b) By L’Hopital’s rule one has 


, t 2 

lim — = 0 

t —>OO g" 

Since /(0) = 0, /(f) is bounded. Since /'(f) = (2f — f 2 )e - *, /(f) has a 
maximum when f = 2. The value of this maximum is /(2) = 4e~ 2 . Hence, 
M = 4e -2 and a — 0 ■ 

Problem 43.20 

Consider the function /(f) = J t+1 . 

(a) Is /(f) continuous on 0 < f < oo, discontinuous but piecewise continuous 
on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 

Solution. 

t 2 

(a) Since ef 2 and e 2 * + l are continuous everywhere, /(f) = J t+1 is continuous 
on 0 < f < oo. 

(b) Since e 2t + 1 < e 2t + e 2t = 2e 2t , f(t) > ^e^e~ 2t = ^e^~ 2t . But for f > 4 

we have f 2 — 2f > C. Hence, /(f) > \e^. So /(f) is not of exponential order 
at infinity ■ 
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Problem 43.21 

Consider the floor function f(t) = |/J> where for any integer n we have 
[_tj = n for all n < t < n + 1. 

(a) Is /(f) continuous on 0 < f < oo, discontinuous but piecewise continuous 
on 0 < f < oo, or neither? 

(b) Are there fixed numbers a and M such that |/(f)| < Me at for 0 < f < oo? 


Solution. 

(a) The floor function is a piecewise continuous function on 0 < f < oo. 

(b) Since |_tj < t < e f for 0 > t < oo we find M = 1 and a — 1 ■ 


Problem 43.22 

Find £-‘ (;y . 


Solution. 

Since C ( —) = —, s > a we find 

\s—a / s—a 1 


c~ 


= 3£ 


-i 


= 3e 2 / t > 0 


Problem 43.23 

Find £-■ (-J + jij) . 

Solution. 

Since C\t I = 4, s > 0 and C (—) = —, s > a we find 

l -i s z 1 \ s—a / s—a 1 



= - 2f + e~\ f > 0 ■ 


Problem 43.24 

Fincl (£ + £)■ 


Solution. 

We have 


£ 


-i 


+ 


s + 2 s — 2 


= 2£ _1 


s + 2 


+ 2£ _1 


s — 2 


= 2(e^ 2t +e 2t ), t > 0 
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Section 44 


Problem 44.1 

Use Table £ to find £[ 2e t + 5]. 

Solution. 

£[2e* + 5] = 2 C[e l ] + 5£[1] = s > 1 ■ 

Problem 44.2 

Use Table £ to find £[e 3t ~ 3 h(t — 1)]. 

Solution. 

C[e 3t ~ 3 h(t — 1)] = C[e 3(t ~^h(t — 1)] = e _s £[e 3< ] = ——s > 3 ■ 

s o 

Problem 44.3 

Use Table £ to find £ [sin 2 cut]. 

Solution. 

£[sin 2 u>t] = £[-— co * 2ujt ] = ^(£[1]—£[cos 2u>t]) = ^ 0 - ^ + ^ 2 ) > s > °" 

Problem 44.4 

Use Table £ to find £[sin3t cos3t]. 

Solution. 

r n r sin 6fl 1 „ r , 3 

£ sm3tcos3f = £ —-— = -£ sm6f = ———, s > 0 ■ 

2 J 2 s 2 + 36 

Problem 44.5 

Use Table £ to find £[e 2t cos3f]. 

Solution. 

£[e a c OS 3i]= (s _ S ~ )2 2 +9 , s > 2 1 
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Problem 44.6 

Use Table C to find C[e 4t (t 2 + 31 + 5)]. 

Solution. 


£[e 4t (l 2 +31+5)] = £[e 4 U 2 ]+3£[e 4i l]+5£[e 4t ] 


2 3 5 
(s-4 ) 3 + (s-4 ) 2 + ^4’ S> " 


Problem 44.7 

Use Table C to find £~ 1 [ s3 ^° 95 + ^]. 

Solution. 


£ 


-ir 


10 


L s 2 + 25 s - 3 


= 2 C 


-ir 


s 2 + 25 


;+4£ 


-ir 


2sin51 + 4e 3t , 1> 0i 


Problem 44.8 

Use Table C to find £~ 1 [ ^_ 5 3 ^ )4 ]. 

Solution. 


C- 1 


5 


(s-3) 



(*- 3) 1 


Uv, 

o 


t > 0 ■ 


Problem 44.9 

Use Table C to find £ _1 [^—4-]. 


Solution. 


C- 4 [^~] = e 9 ^h(t - 2 ) 

s — 9 


0 , 0 < 1 < 2 
e 9 ( t_2 ), 1 > 2 ■ 


Problem 44.10 

Use Table C to find C~ l [ e 


Solution. 

We have 


C 


-i 


e“ 3s (2s + 7)- 
s 2 + 16 


=2 C 


-H 


e 3s s 


+ U- 1 ' 


-3*4 


s 2 + 16 4 s 2 + 16 


7 

=2 cos 4(1 — 3 )h(t — 3) + - sin 4(1 


- 3 )h(t — 3), 1 > 0 
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Problem 44.11 

Graph the function /(f) = h(t — 1) + h(t — 3) for f > 0, where h(t) is the 
Heaviside step function, and use Table £ to find £[/(£)]. 

Solution. 

Note that 

( 0, 0 < t < 1 
f(t) = { 1, 1 < t < 3 
[2, t > 3 

The graph of f(t ) is shown below. Using Table £ we find 

_ g _ 05 

U[/(t)] = £[h(t - 1)] + C[h(t - 3)] = + e —^~ , s > 0 ■ 


2 




1 


3 


* 


Problem 44.12 

Graph the function f(t ) = f[/i(t — 1) — /r(t — 3)] for t > 0, where h(t ) is the 
Heaviside step function, and use Table C to find £[/(£)]. 

Solution. 

Note that 

f 0, 0 < t < 1 
/(t) = ^ t, 1 < t < 3 
[0, t > 3 

The graph of /(f) is shown below. Using Table £ we find 

£[/(t)] =£[(£ — l)h(f — 1) + h(t — 1) — (f — 3 )h(t — 3) — 3h(f — 3)] 

=£[(f — l)h(f — 1)] + £[h(t — 1)] — £[(t — 3 )h(t — 3)] — 3 £[h(t — 3)] 

e ~s e -s e -3s 3 e - 3s 

= —tt H-n-, s > 1 ■ 
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3 
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Problem 44.13 

Graph the function /(f) = 3 [h(t — 1) — h(t — 4)] for t > 0, where h(t) is the 
Heaviside step function, and use Table £ to find £[/(f)]. 


Solution. 

Note that 

( 0, 0 < f < 1 
fit) = l 3, 1 < f < 4 
[ 0, f > 4 

The graph of /(f) is shown below. Using Table C we find 

Q P -s Se _4s 

£[/(f)] = 3 C[h{t — 1)] — 3 jC[h(t — 4)1 =-, s > 0 ■ 

s s 



1 



*► 


Problem 44.14 

Graph the function /(f) — \2 — t\[h(t — 1) — h(t — 3)] for f > 0, where h(t) is 
the Heaviside step function, and use Table C to find £[/(f)]. 

Solution. 

Note that 

( 0, 0 < f < 1 

fit) = l 12 — f|, 1 < f < 3 
I 0, t > 3 
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The graph of /(£) is shown below. Using Table C we find 


£[/(£)] =(2 - £)/i(f - 1) + 2(f - 2 )h(t - 2) - (t - 2)h(t - 3) 

=£[-(£ - 1 )h(t - 1) + h(t - 1) + 2(£ - 2 )h(t - 2) - (£ - 3 )h(t - 3) - h(t - 3)] 
= - £[(t - 1 )h(t - 1)] + C[h{t - 1)] + 2 C[(t - 2 )h(t - 2)] 

-£[{t - 3 )h(t - 3)] - £[h(t - 3)] 

e~ s e~ s 2e~ 2s e" 3s e" 3s 
— — — 7T T - T ^— — 7T~ — -, S > 0 | 


1 





Problem 44.15 

Graph the function f(t) = h(2 — t) for t > 0, where h(t) is the Heaviside step 
function, and use Table C to find £[f(t)]. 


Solution. 

Note that 


m = 


1, 0 < t < 2 
0, t > 2 


The graph of /(f) is shown below. From this graph we see that /(f) = 
h(t) — h(t — 2 )h{t — 2). Using Table C we find 


£[/(£)] = £[h(t)} - £[h(t - 1 )h(t - 1)] 



s > 0 ■ 
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Problem 44.16 

Graph the function /(f) = h(t — 1) + h( 4 — f) for t > 0, where h(t) is the 
Heaviside step function, and use Table £ to find £[/(f)]. 

Solution. 

Note that 

( 1, 0 < t < 1 
fit) = \ 2, 1 < t < 4 
( 1, ^>4 

The graph of /(f) is shown below. Using Table £ we find 

—s /*4 I | —s —4s 

£[/(f)] = £[h{t— l)]+£[/i(4—f)] = --b / e~ st dt =---, s > 0i 

S ,/n S 


2 




4 


Problem 44.17 

The graph of /(f) is given below. Represent /(f) as a combination of Heav¬ 
iside step functions, and use Table £ to calculate the Laplace transform of 

fit)- 



Solution. 

From the graph we see that 

fit) = (f — 2 )h(t — 2) — (f — 3)/i(f — 3) — /i(f — 4) 

Thus, 

_ 2 § _3s _4s 

£[/(f)] = £[(f—2)/i(f—2)1—£[(f—3)/i(f—3)1—£[/i(f—4)1 =- — --, s > 0■ 

s 2 s 
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Problem 44.18 

The graph of f(t ) is given below. Represent f(t ) as a combination of Heav¬ 
iside step functions, and use Table C to calculate the Laplace transform of 
/(*)■ 


0 1 


2 



Solution. 

From the graph we see that 

fit) = hit — 1) + hit — 2) — 2 h(t — 3). 


Thus, 

C[f(t)] = C[h(t-l)]-2C[h(t-3)] + C[h(t-2)] = - 2( ' ( 2 \ s > 0 ■ 

Problem 44.19 

Using the partial fraction decomposition find £ -1 ( s _ 3 )^ s+1 ) • 

Solution. 

Write 

12 _ A B 

(s — 3)(s + 1) s — 3 s T 1 

Multiply both sides of this equation by s — 3 and cancel common factors to 
obtain 

12 „ Bis-3) 

-= A H--- 

s+1 s+1 

Now, find A by setting s = 3 to obtain A = 3. Similarly, by multiplying both 
sides by s + 1 and then setting s — — 1 in the resulting equation leads to 
B = —3. Hence, 


12 ( 1 1 ^ 

(s — 3)(s + 1) ys — 3 s + 1 J 
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Finally, 


£- 


12 


L(s - 3)(s + 1) 


=3£- 1 

1 

-3 zr 1 

i 

_s — 3 

+ 1 




=3e M - 3e , t > 0 


Problem 44.20 

Using the partial fraction decomposition find ZU 1 


24e~ 5s 
s 2 —9 


Solution. 

Write 

24 A | B 

(s-3)(s + 3) ~ s — 3 + s + 3" 

Multiply both sides of this equation by s — 3 and cancel common factors to 
obtain 

24 , B(s- 3) 

-= A H- 2 - 

s + 3 s + 3 

Now, find A by setting s = 3 to obtain A = 4. Similarly, by multiplying both 
sides by s + 3 and then setting s = —3 in the resulting equation leads to 
B = —4. Hence, 


_Hi_ = 4 fM_ 

(s-3)(s + 3) V s - 3 


s + 3 


Finally, 


ZT 1 



24e- 5s 
-3)(s + 3) 


= 4£ _1 

- g — 5s ' 

1 

- e -5s - 

_s — 3 

s + 3 




=4[ e 3(t - 5) - e~ z(t - b) }h(t - 5), t> 0 ■ 


Problem 44.21 

Use Laplace transform technique to solve the initial value problem 


where 


y' + 4y = g(t), y( 0) = 2 


| 0, 0 < t < 1 
g{t) =\ 12, 1 < t < 3 
I 0, t > 3 
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Solution. 

Note first that g(t) = 12 [h(t — 1) — h{t — 3)] so that 


C[g(t)] = 12 C[h{t - 1)] - 12 C[h(t - 3)] = 12(e *— 6 s > 0. 
Now taking the Laplace transform of the DE and using linearity we find 

C[y'} +4 C[y] = C[g(t)\. 

But C[y'] = s£[y\ — y( 0) = sC[y\\ — 2. Letting £[y\ = Y(s) we obtain 


sy (s) - 2 + 4Y(s) = 12 


e s — e 3s 


Solving for Y(s) we End 


Y(s) = 


12 - 


e s — e 3s 


But 


and 


£ 


-l 


s + 4 s(s + 4) 
2 


s + 4 


= 2e 


—At 


c 


-1 


12 - 


e s — e 3s 
s(s + 4) 


=3 C 


-i 


(e~ s - e~ 3s ) 


1 


=3C 


-i 


3C 


-l 


s s + 4 

' -3s 


3C 


-i 


s + 4 


+ 3£ 


-l 


3s 


s + 4 


=3 hit — 1) — 3 h(t — 3) — 3e l ^h(t — 1) + 3e ^ 3 ^h(t — 3) 


ffence, 

y(t) = 2e~ u +3[h(t—l)—h(t—3)] — 3[e~ 4 ^ 1 ^h(t—l)—e~ 4: ^~ 3 ^h(t—3)], t > 0 ■ 

Problem 44.22 

Use Laplace transform technique to solve the initial value problem 
y"-4y = e 3 \ 2 /( 0 ) = 0, j/(0) = 0. 
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Solution. 

Taking the Laplace transform of the DE and using linearity we find 


£[y"}-A£[y] = C[e 3t }. 

But C[y"] = s 2 C[y\ — sy(0) — y'( 0) = s 2 C[y). Letting C[y\ = Y (s) we obtain 

s 2 y( s )-4y( s ) = -^. 

S — 3 


Solving for Y (s) we find 

Y(8) = 


(s — 3)(s — 2)(s + 2)' 

Using partial fraction decomposition 

1 ABC 


(s — 3)(s — 2)(s T 2) s — 3 s + 2 s — 2 
we End A = B — A, and C = — Thus, 


y(*) =£ 1 

l 


L(s — 3)(s — 2)(s + 2) 


=w 


s — 3 


-I_ 

+ 20 


s + 2 


= 2 e 3i i JL e -2t _ I e » t> n 

5 + 20 4 e ’ - U 



1 

s-2 


Problem 44.23 

Obtain the Laplace transform of the function f* f(X)dX in terms of £[/(£)] = 
F(s) given that / Q 2 f(X)dX = 3. 


Solution. 

We have 


C 


f(X)dX 


IJ 2 


=£ 




/(A)dA - / f(\)d\ 
Jo 

-m 


F(s) 3 
■ 7 s > 0 


s s 



Section 45 


In Problems 45.1 - 45.4, give the form of the partial fraction expansion for 
F(s). You need not evaluate the constants in the expansion. However, if the 
denominator has an irreducible quadratic expression then use the completing 
the square process to write it as the sum/difference of two squares. 

Problem 45.1 


F(s) = 

(s — l) 3 (s — 2) 2 ’ 

Solution. 

T-i/ \ _ -4-i -4-2 ^3 Bi B-2 

F ^-T^Ty + (^Ty + — + T^F + —. a 

Problem 45.2 


Solution. 


Problem 45.3 


F(s) 


s 2 + 5s - 3 
(s 2 + 16) (s — 2) ’ 


F(s) 


A\s + A 2 B\ 
s 2 + 16 + s-2 " 


F(s) = _-_-_ 

{) (s 2 + l) 2 (s + 4) 2 ’ 

Solution. 

= + A 2 , ^3-S + A 4 B x B 2 

1J (s 2 + l) 2 s 2 + 1 (s + 4) 2 s + 4 ■ 

Problem 45.4 


s 4 + 5s 2 + 2s - 9 
( ' S) ~ (s 2 + 8s + 17)(s — 2) 2 ’ 
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Solution. 


t-,/ \ _ Ai A 2 B\s + B'2 

( ' S) ~ (s - 2) 2 + ^2 + (s + 4) 2 + 1 " 


Problem_45.5 

Find C~ 


(s+l ) 3 


Solution. 

Using Tabic C we find £ _1 


1 

(s+ 1) 3 


|e B 2 , t > 0 ■ 


Problem 45.6 

Find [;£&] ■ 

Solution. 

We factor the denominator and split the rational function into partial frac¬ 
tions: 


2s - 3 A | B 

(s - l)(s - 2) “ + 1^2’ 

Multiplying both sides by (s — l)(s — 2) and simplifying to obtain 

2s - 3 =ZL(s - 2) + B(s - 1) 

= (A + B)s -2 A- B. 


Equating coefficients of like powers of s we obtain the system 

f A+B = 2 

I -2A-B = -3. 


Solving this system by elimination we find A = 1 and B — 1. Now finding 
the inverse Laplace transform to obtain 


£ 


-1 


2s-3 


L(s-1)(s-2) 


= zr 1 

1 

+ C^ 

1 

s — 1 

_s — 2_ 




= e l + e 2t , t > 0 


Problem 45.7 


Find C~ 


4s 2 +S+l 

s 3 +s 



Solution. 

We factor the denominator and split the rational function into partial frac¬ 
tions: 


4s 2 + s + l A Bs A C 
s(s 2 + 1) s + s 2 + 1 

Multiplying both sides by s(s 2 + 1) and simplifying to obtain 

4s 2 + s A 1 =A(s 2 + 1) + (Bs A C)s 
=(A + B)s 2 + Cs + A. 


Equating coefficients of like powers of s we obtain A A B = 4, C — 1, A = 1. 
Thus, B — 3. Now finding the inverse Laplace transform to obtain 

x I” 4s 2 + s + 1 
[ s(s 2 + l) 

=1 + 3 cos t A sin t, t > 0 ■ 


=C 


-l 


+ 3 C~ l 


s 2 + 1 


AC 


-i 


s 2 + 1 


Problem 45.8 


Find C~ 


++6s+8 

s 4 +8s 2 +16 


Solution. 

We factor the denominator and split the rational function into partial frac¬ 
tions: 


s 2 A 6s + 8 B\s A Ci B 2 s A C 2 

(s 2 A 4) 2 s 2 A 4 (s 2 + 4) 2 

Multiplying both sides by (s 2 + 4) 2 and simplifying to obtain 

s" A 6s + 8 =(B\S + Ci)(s 2 + 4) + B 2 s + C 2 

=f?is 3 + Li'S 2 + (4-Bi + B 2 )s A 4 Ci A C 2 . 


Equating coefficients of like powers of s we obtain Bi = 0, C\ — 1 ,B 2 = 6, 
and C 2 = 4. Now finding the inverse Laplace transform to obtain 


C l 


s 2 A 6s + 8 
(s 2 + 4) 2 


=£ 


-i 


s 2 A 4 


+ 6 C 


-l 


L(s 2 + 4)^ 


+ 4 C 


-i 


(s 2 + 4)^ 


= - sin 2t A 6 ( - sin 2t ) + 4 ( — [sin 2 1 — 2 1 cos 2 1] 

2 \4 / \16 

3 3 1 

= -t sin 2 1 A - sin 2 1 - 1 cos 2 1, t > 0 ■ 

2 4 2 



Problem 45.9 

Use Laplace transform to solve the initial value problem 

y + 2y = 26 sin 3 1, y( 0) = 3. 


Solution. 

Taking the Laplace of both sides to obtain 

C[y] + 2 C[y\ = 26£[sin 3f]. 

Using Table C the last equation reduces to 

sY(s) - y(0) + 2 Y(s) = 26 ■ 

Solving this equation for Y(s) we hnd 

V7U - 3 I 78 

“ s + 2 + (s + 2)(s 2 + 9)' 

Using the partial fraction decomposition we can write 

1 _ A Bs + C 

(s + 2)(s 2 + 9) “ s + 2 + s 2 + 9 ‘ 

Multipliying both sides by (s + 2)(s 2 + 9) to obtain 

1 =A(s 2 + 9) + (yBs + C)(s + 2) 

—{A + B)s 2 + (2 B + C)s + 9A + 2 C. 


Equating coefficients of like powers of s we hnd A + B = 0, 2 B + C = 0, and 
9 A + 2U = 1. Solving this system we hnd A = - jb, B = — jb, and C — A. 
Thus, 

y(s) = db - 6 (^Tq) + 4 (iJTg)' 

Finally, 


y(t) = = 9/r 1 

= 9e 


— 6 £ 


-i 


_s + 2 

~ 2< — 6 cos 3t + 4 sin 3 t, t > 0 


s 2 + 9 


+ 4£ 


-i 


s 2 + 9 
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Problem 45.10 

Use Laplace transform to solve the initial value problem 

y' + 2y = 4 1, y{ 0) = 3. 


Solution. 

Taking the Laplace of both sides to obtain 

C[y'] + 2C[y\ = 4 £[t]. 

Using Table £ the last equation reduces to 

5 y(s)-j/(o) + 2 y( s ) = 4- 

s z 

Solving this equation for Y(s) we hnd 

U s + 2 (s + 2)s 2 ' 

Using the partial fraction decomposition we can write 

1 _ A Bs + C 

(s + 2)s 2 s + 2 + s 2 

Multipliying both sides by (s + 2)s 2 to obtain 

1 = As 2 + (Us + C)(s + 2) 

= (A + B)s 2 + (2 B + C)s + 2 C 


Equating coefficients of like powers of s we hnd A + B = 0, 2 B + C = 0, and 
2U = 1. Solving this system we hnd A — B = — and C — |. Thus, 


Y{s) 


4 1 1 

s + 2 s + s 2 


Finally, 


y(t) 


jC-ifyfs)] = 4C~ X — 

Ls + 

4e~ 2t -1 + 2 1, t> 0 ■ 



+ 2£ _1 
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Problem 45.11 

Use Laplace transform to solve the initial value problem 


y" + 3 y' + 2 y = 6e~% y{ 0) = 1, j/'(0) = 2. 

Solution. 

Taking the Laplace of both sides to obtain 

£{y"}+3£{y'}+2£[y} = 6£[e- t }. 

Using Table £ the last equation reduces to 

s 2 Y(s) - sy(0) - y'(0) + 3(sU(s) - 2/(0)) + 2Y(s) = 

s + 1 

Solving this equation for Y(s) we hnd 

s + 5 6 s 2 + 6s + 11 

(S) = (s + l)(s + 2) + (s + 2)(s + l) 2 = (s + l) 2 (s + 2) ’ 

Using the partial fraction decomposition we can write 

s 2 T 6s T11 A B C 

(s + 2)(s + l) 2 “ s + 2 + s + 1 + (s + l) 2 ' 

Multipliying both sides by (s + 2)(s + l) 2 to obtain 

s 2 + 6s + 11 = A(s + l) 2 + B(s + l)(s + 2) + C(s + 2) 

= (A + B)s 2 + (2 A + 3 B + C)s + A + 2B + 2C 

Equating coefficients of like powers of s we hnd A + B = 1, 2A + ?>B + C = 6, 

and A + 2 B + 2 C = 11. Solving this system we hnd A = 3, B — —2, and 

C — 6. Thus, 

v7-/\ 3 2 6 

( ' S) “ s + 2 s + 1 + (s + l) 2 ‘ 

Finally, 

y(t) = £~ 1 [Y(s)] = 3£- 1 —— - 2/U 1 —— + 6/T 1 - 1 
W L V ;J [s + 2j [s + lj [(s + l) 2 J 

= 3e^ 2 * - 2e~ f + 6 te^, t > 0 ■ 
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Problem 45.12 

Use Laplace transform to solve the initial value problem 
lj" + % = cos 2 1, y{ 0) = 1, y'{ 0) = 1. 


Solution. 

Taking the Laplace of both sides to obtain 

C[y"] + 4 C[y\ = £[cos2t]. 

Using Table C the last equation reduces to 

s 2 U(s) - sy( 0) - ?/(0) + 4 Y(s) = 

s z + 4 


Solving this equation for Y(s) we find 

s + 1 


Y(*) = 


s 2 + 4 (s 2 + 4) 2 


Using Table C again we find 

s 


y(t) = £ 


-1 


s 2 + 4 


I ^- 1 _ 

2 s 2 + 4 


ZU 1 


L (, 2 


1 t 

= cos 2t + - sin 2t + - sin 2 t, t > 0 


Problem 45.13 

Lise Laplace transform to solve the initial value problem 
y"-2y' + y = e 2t , y( 0) = 0, j/'(0) = 0. 


Solution. 

Taking the Laplace of both sides to obtain 

C[y"]-2C[y']+C[y\=C[e 2t \. 
Using Table C the last equation reduces to 


s 2 U(s) - S y(0) - y'{ 0) - 2sU(s) + 2y(0) + y(s) 


s 

+ 4p 


1 

s — 2' 
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Solving this equation for Y(s) we find 


Y{s) = 


(■s - 1) 2 (® - 2) ■ 

Using the partial fraction decomposition, we can write 

Y(a) = * B C 


s-1 (s-1) 2 s — 2' 

Multipliying both sides by (s — 2)(s — l) 2 to obtain 

1 = A(s - l)(s - 2) + B(s - 2) + C(s - l) 2 
= (A + C)s 2 + (-3 A + B - 2C)s + 2A-2B + C 

Equating coefficients of like powers of s we find A+C = 0, —3A+B — 2C = 0, 
and 2 A — 2 B + C — 1. Solving this system we find A = — 1, B — — 1, and 
C — 1. Thus, 

I'M = "A - 7 1 1 


Finally, 

r/(t) = U _1 [U(s)] = -CT 1 


l) s 




l) 2 


+ U" 1 


= —e 


te* + e 2t , t > 0 


Problem 45.14 

Use Laplace transform to solve the initial value problem 

y" + 9y = g(t), y( 0) = l, j/(o) = 3 

where 

a(t) = { * 


0 < t < 71 
71 < t < OO 


Solution. 

Taking the Laplace of both sides to obtain 

C[y"] + 9 C[y\ = C[g(t )] = 6 C[h{t) - h(t - 7r)]. 
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Using Table C the last equation reduces to 


s 2 l"(s) - sj/(0) - p'(0) + 9U(s) = - 


6 6e 7rs 


Solving this equation for Y(s) we hnd 


Y(s) = 


s + 3 


6 


-(1 — e _7rs ). 


s 2 + 9 s(s 2 + 9) 

Using the partial fraction decomposition, we can write 

6 A Bs + C 
s(s 2 + 9) s + s 2 + 9 

Multiplying both sides by s(s 2 + 9) to obtain 


6 = A(s 2 + 9) + (Bs + C)s 
= (A + B)s 2 + Cs + 9A 


Equating coefficients of like powers of s we hnd A + B = 0, C — 0, and 
9 A = 6. Solving this system we hnd A — B = — and C — 0. Thus, 


Y{s) 


s 2 + 9 


+ 


3 

s 2 + 9 


+ (1 - e~ ns ) 



2 s \ 

3 s 2 + 9 ) ' 


Finally, 

2 2 

y(t) = £ _1 [y'( s )] = cos 3 1 + sin3f + -(1 — cos 3 1) — -(1 — cos 3 (t — n))h(t — it) 

o o 

2 2 

= cos 3 1 + sin 3 1 + -(1 — cos 3 1 ) — -(1 + cos 3 t)h(t — vr), t > 0 ■ 
o o 


Problem 45.15 

Determine the constants a, (3, tjo, and y' 0 so that Y (s) = is the Laplace 

transform of the solution to the initial value problem 


y" + ocj/ + f3y = 0, 2/(0) = 2 / 0 , 2/(0) = y'o- 


Solution. 

Taking the Laplace transform of both sides we hnd 

s 2 U(s) - s2/o - y' 0 + asY (s) - ay G + f$Y(s) = 0. 
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Solving for Y(s) we find 


Y( \ = s Vo + (y'o + a Vo) = 2s - 1 
s 2 + as + (3 s 2 + s + 2 

By identification we find a — 1, (3 = 2, y 0 = 2, and t/q = — 3 ■ 

Problem 45.16 

Determine the constants a,/3,yo, and i/ 0 so that Y(s) = is the Laplace 
transform of the solution to the initial value problem 

y" + ay' + (3y = 0, y{ 0) = y 0 , y\ 0) = y' 0 - 


Solution. 

Taking the Laplace transform of both sides we find 

s 2 Y(s) - sy 0 -y' 0 + asY(s) - ay 0 + (3Y(s) = 0. 

Solving for Y(s) we find 

A = sy 0 + (y' 0 + ay 0 ) = s 

s 2 + as + (3 s 2 T 2s T 1' 

By identification we find a = 2, (3 = 1, yo = 1, and t/q = — 2 ■ 
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Section 46 

Problem 46.1 

Find the Laplace transform of the periodic function whose graph is shown. 

-0 


12 3 4 



Solution. 

The function is of period T = 2. Thus, 


_ — st 


dt + 


— St 


dt 



Hence, 

3 - 2e~ s - e~ 2s 
s(l — e ~ 2s ) * 



-(3 — 2e~ s — e~ 2s ). 
s 


Problem 46.2 

Find the Laplace transform of the periodic function whose graph is shown. 



Solution. 

The function is of period T = 4. Thus, 


—st 


dt + 


—st 


dt 



— st 


-I 3 


= -(2 - e~ s - e~ 6s ). 
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Hence, 


£[/(*)] 


2 — 6 S — 6 
s(l — e ~ 4s ) * 


Problem 46.3 

Find the Laplace transform of the periodic function whose graph is shown. 



Solution. 

The function is of period T = 2. Thus, 


(t — l)e st dt — I te st dt — e st dt 


t _ f e st e st 

-—e - - —I- 

s s z s 

= -?[(l + ^)e- s -l] 


Hence, 


cm] = 


s 2 (l — e _2s ) 


[i_ (s+ i )e -.] 


Problem 46.4 

Find the Laplace transform of the periodic function whose graph is shown. 
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Solution. 

The function is of period T = 2. Thus, 


te st dt = 


~ — {st + l)e 


—st 


= --[(2s+ 1 )e- 2 *- 1 ] 


Hence, 


mt)] = 


s 2 (l — e _2s ) 


[1 - (2s + l)e~ 2s ) 


Problem 46.5 

State the period of the function f{t) and find its Laplace transform where 


m = 


sin t, 0 < t < 7T 
0, 7T < t < 2 tt 


/(t + 2vr) = fit), t > 0. 


Solution. 

The graph of f(t) is shown below. 



The function f(t) is of period T = 2n. The Laplace transform of f(t) is 

sin te~ st dt 
1 - e _2?rs 



Using integration by parts twice we find 

/* iD — St 

-St, 


sinte st dt = 


1 + s 2 


(cos t + s sin t) 
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Thus, 


/ 

Jo 


sin te st dt = -- (cos t + s sin t ) 

1 + s z 

e-ns x 


0 


1 + s 2 1 + s 2 
1 + e~ ns 
1 + s 2 


Hence, 


(1 + s 2 )(l — e~ 2ns ) * 


1 + 


Problem 46.6 

State the period of the function f(t) = 1 — e _t , 0 < t < 2, f(t + 2) = f(t), 
and find its Laplace transform. 

Solution. 

The graph of f(t) is shown below 



2 4 


The function is periodic of period T = 2. Its Laplace transform is 



But 


/ 


(1 -e-^e-^dt = 


_st ] 2 r e -(s+ 1 )*i 2 2 



s +1 


Hence, 
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Problem 46.7 

Using Example 46.3 find 


C~ l 


s 2 — s e s 

s 3 + s( 1 - e ~ s ) 


Solution. 

Note first that 

s 2 - s e~ s _ 1 /I se~ s \ 

s 3 + s(l-e- s ) s \s 2 s 2 (l — e~ s ) J 

Using Example 23.3, we find 

f(t) = 1 -git) 

where g{t) is the sawtooth function shown below 


m 



Problem 46.8 

An object having mass m is initially at rest on a frictionless horizontal surface. 
At time t = 0, a periodic force is applied horizontally to the object, causing 
it to move in the positive x-direction. The force, in newtons, is given by 

r f 0 , o<t<f 

m = { f(t+r) = f(t), t> o. 

[ 0, I < t < T 


101 












The initial value problem for the horizontal position, x(t), of the object is 

mx"(t ) = /(f), x(0) = a/(0) = 0. 

(a) Use Laplace transforms to determine the velocity, v(t) = x'(t), and the 
position, x(t), of the object. 

(b) Let m — 1 kg, f 0 = 1 N, and T = 1 sec. What is the velocity, v, and 
position, x, of the object at t — 1.25 sec? 


Solution. 

(a) Taking Laplace transform of both sides we find ms 2 X (s) = /0 'j 1 / 
. Solving for X (s) we find 

1 


_ fo e Btdt 

sT — 


fo ( l-e~ s ^ 
s l—e~ sT 


Y( \ f° 

A ( s ) = 


m s° i -f e - s l 


Also, 


V(s) = £[«(()] = sX(s) = k ■ 4 ■ 1 


1 1 


m s 4 l + e" s ? m s S (l + e - s t) 
Hence, by Example 46.1 and Table £ we can write 


v(t) = — / f(u)du. 


m 


Since X(s) = 


= * f( f )} we have 

x(t) = — (f * /(f)) = — f (f - u)f(u)du. 
m m In 


(b) We have x(1.25) = ,/ 0 2 (f — u)du + f 1 4 (^—w)dw = A meters and u(1.25) = 

5 15 

f 0 4 fi u )du = f 0 2 dt + J 4 dt = | m/sec ■ 


Problem 46.9 

Consider the initial value problem 


ay" + by' + cy = /(f), 7/(0) = t/( 0) = 0, f > 0. 

Suppose that the transfer function of this system is given by <f>(s) = 9 5 2 + 1 5s+ 9 

(a) What are the constants a, b , and c? 

(b) If /(f) = e~ l , determine F(s), Y(s), and y(t). 
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Solution. 

(a) Taking the Laplace transform of both sides we find as 2 Y(s) + bsY(s) + 
cy(s) = F(s) or 


$(s) 


Y{s) 

F(s) 


1 _ 1 
as 2 + bs + c 2s 2 + 5s + 2 


By identification we hnd a = 2, 6 = 5, and c = 2. 
(b) If f{t) = e _t then F(s) = £[e _t ] = -A-p Thus, 


y(s) = 4>(s)F(s) 


1 

(s + l)(2s 2 + 5s + 2)' 


Using partial fraction decomposition 

1 _ A , B < C 
(s + l)(2s + l)(s + 2) ““ s + 1 + 2s + 1 + s + 2 

Multiplying both sides by s + 1 and setting s = — 1 we find A = — 1. Next, 
multiplying both sides by 2s+ 1 and setting s = — | we find B = Similarly, 
multiplying both sides by s + 2 and setting s = —2 we hnd C — Thus, 


3/(f) = - £ 1 


1 

+ \ c ~ l 

1 

+ 

1 

_s + 1 

- s + 

s + 2 

3 

3 


_+ * _ 
= — e + -e 
3 


5 + -e" 2t ,f > 0 
o 


Problem 46.10 

Consider the initial value problem 


ay" + by' + cy= /(f), y(0) = y'(0) = 0, t > 0 


Suppose that an input /(f) = f, when applied to the above system produces 
the output y(t) = 2(e _t — 1) + f(e+ 1), t > 0. 

(a) What is the system transfer function? 

(b) What will be the output if the Heaviside unit step function /(f) = h(t) 
is applied to the system? 


Solution. 

(a) Since /(f) = f we hnd F(s) = Aslo, y(s) = C[y(t)] = £[2e _ * — 2 + 
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te“‘+i] = A 

If /(f) = hit) 


~s + 


+ 4 = 


2 (s+l ) 2 


(s+1) 2 

(b) If /(f) = h(t) then F(s) = ^ and Y(s) = 
partial fraction decomposition we find 


But 4ms) = 


_ Y(s) _ 1 


F(s) 

4>(s)F(s) 


ww 

1 

s(s+l ) 2 


Using 


1 _A B C 
s(s + l) 2 s s + 1 (s + l) 2 

1 =A(s + l) 2 + Bs(s + 1) + Us 
1 =(A + B)s 2 + (2 A + B + C)s + A 


Equating coefficients of like powers of s we find A = 1, B = — 1, and 
U = — 1. Therefore, 


and 


Y{s) = 



1 

s + 1 


1 

(s + 1) 2 


y(t) = C 1 |Y(s)] = 1 — e 4 — te *, f > 0 ■ 


Problem 46.11 

Consider the initial value problem 

y" + j/ + y = m, y(0)=y'(0) = 0, 


where 


f 1, 0 < f < 1 

f(t) = { f(t + 2) = f{t ) 

1-1, 1 < f < 2 


(a) Determine the system transfer function 4>(s). 

(b) Determine Y(s). 


Solution. 

(a) Taking the Laplace transform of both sides we hnd 

s 2 Y (s) + sY (s) + Y (s) = F(s) 


so that 


4>(s) = 


Y(s) 

F(s) 


1 

^2 j g j ^ 
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(b) But 


Hence, 


and 


f(t)e st dt 




e st dt 


" ^ — st - 

1 

" ^ — st " 

e 


e 

—s 

0 

—s 


= -(l-e~ s ) + -(e- 2s -e~ s ) 
s s 

_(1 -e“ s ) 2 


s 


F(s) 


( 1 - e -*) 2 
s(l — e~ 2s ) 


(i - g—) 
s(l + e~ s ) 


Y(s) = $(s)F(s) 


(!-«-) _ 

s(l + e -s )(s 2 + s + 1) 


Problem 46.12 

Consider the initial value problem 


y'" - % = e f + t f 2/(0) = 2/(0) = 2/"(0) = 0. 


(a) Determine the system transfer function <l>(s). 

(b) Determine Y(s). 


Solution. 

(a) Taking Laplace transform of both sides we find 

s 3 Y(s) -4Y(s) = F(s). 


Thus, 


(b) We have 


Hence, 


4>(s) 


Y{s) 
F(s) 


1 


s 3 — 4 


F(s) - C[e f + t] - j—j + ^ 


s 2 + s — 1 
(s - l)s 2 ■ 


Y{s) 


s 2 + s — 1 
s 2 (s — l)(s 3 — 4) * 
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Problem 46.13 

Consider the initial value problem 

y" + W + cy= 6(f), y( 0) = y 0 , y'( 0) = y' 0 , t> 0. 

Suppose that C[y(t)\ = Y(s) = . Determine the constants 6, c, yo, 

and y' 0 . 

Solution. 

Take the Laplace transform of both sides to obtain 

X 

s 2 F(s) - sy 0 ~y' 0 + bsY(s) - by 0 + cY'(s) = 

Solve to find 


Y( , _ s 2 yo + s(y' 0 + by 0 ) + 1 
s 3 + bs 2 + cs 
s 2 + 2s + 1 
~s 3 + 3s 2 + 2s' 

By comparison we find 6 = 3, c = 2, yo — 1, and y' 0 + byo = 2 or y' 0 = — 1 ■ 
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Section 47 


Problem 47.1 

Consider the functions /(f) — g(t) — h(t), t > 0 where hit) is the Heaviside 
unit step function. Compute f * g in two different ways. 

(a) By directly evaluating the integral. 

(b) By computing £~ 1 [F(s)G(s)] where F(s) = £[/(t)] and G(s) = C[g(t)]. 

Solution. 

(a) We have 

(f*9)(t)= I f(t-s)g(s)ds = 

Jo 

(b) Since F(s) = G{s ) = C[h(t)] = ^ we have (/ * g)(t) = £ _1 [-F(s)Cr(s)] = 

C~ 1 [j i ]=t,t>Qu 

Problem 47.2 

Consider the functions /(f) = e* and g{t) = e~ 2t , t > 0. Compute f * g in 
two different ways. 

(a) By directly evaluating the integral. 

(b) By computing £ _1 [F(s)G(s)] where F(s) = £[/(t)] and G(s) = C[g(t)]. 


h(t — s)h(s)ds = / ds — t, t > 0. 


Solution. 

(a) We have 


(f*9)(t)= / f(t ~ s)g(s)ds = / e (t s) e 2s ds 
Jo Jo 


—e 


e~ 3s ds = 


e* — g 


,t> 0. 

o 

(b) Since F(s) = £[e/ = ^ and G(s) = C[e ~ 2t ] = ^ we find (/ * g)(t) = 
£ _1 [F(s)G(s)] = £~ 1 [ ^ s _ 1 j 1 (s _ 2) ]. Using partial fractions decomposition we 
find 

- - -= -(—-—). 

(s — l)(s + 2 ) 3 v s — 1 s + 2 
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Thus, 


1 / i i \ P t _ P -2t 

(/*<?)(*) = £-^{ 8 ) 0 ( 8 )] = - ^- 1 [—] - /T 1 )—]J = - 3 -,« > 01 

Problem 47.3 

Consider the functions /(f) = sinf and g{t ) = cost, f > 0. Compute f * g in 
two different ways. 

(a) By directly evaluating the integral. 

(b) By computing £ -1 [.F(s)G(s)] where F(s) = £[/(£)] and G(s) = C[g(t)}. 

Solution. 

(a) Using the trigonometric identity 2 sin p cos q = sin (p + q) +sin (p — q ) we 
find that 2 sin (f — s) cos s = sin f + sin (f — 2s). Hence, 

(f * g)(t) = / f(t — s)g(s)ds — / sin (t — s) cos sds 

Jo Jo 

= -[ / sin tds + / sin (f — 2s)ds] 

2 io Jo 

tsinf If 1 . 

=-(— / sin wait 

2 4 J_ t 

t sin f 

(b) Since F(s) = £[sinf] = -^yy and G(s) = £[cosf] = we find 

( f*9)(t ) = £- 1 [F(s)G(s)] = £~ x [ + ^ 2 ] = ^ sinf, t > 0 ■ 

Problem 47.4 

Compute and graph f * g where /(f) = h(t) and g(t) = t[h(t) — h(t — 2)]. 

Solution. 

Since /(f) = h(t), F(s ) = U Similarly, since g{t) = ffi(f) — th(t — 2) = 
ffi(f) - (f - 2)h(t - 2) - 2h{t - 2), G(s) = £ - ^ Thus, F(s)G(s) = 

— ^~s~ — 2e J ■ It follows that 

S ' 3 s° s z 

(/ » g)(t) = ^ - 2) - 2(C - 2)A(« - 2),i > 0. 

The graph of (/ * g)(f) is given below ■ 
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Problem 47.5 

Compute and graph f * g where f(t) = h(t ) — h(t — 1) and g(t) = h(t — 1) 
2h(t — 2). 


Solution. 

Since f(t ) = h(t) — h(t — 1), F(s ) = ^ — Similarly, since g(t) = h(t 
1) - 2 h(t - 2), G(s) = ^ Thus, 


It follows that 


F( S )G( S ) 


e s - 3e 2s + 2e 3s 



(/ * g )(t) = (t - 1 )h(t - 1) - 3 {t - 2 )h(t - 2) + 2(t - 3 )h(t -3 ),t> 0. 
The graph of (/ * g)(t ) is given below ■ 
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Problem 47.6 

Compute t*t*t. 


Solution. 

By the convolution theorem we have £[t*t*t] = (£[/:]) 3 = Hence, 

t*t*t = £~ 1 [j s ] = Ti = m>t>0u 


Problem 47.7 

Compute h(t) * e~ f * e~ 2t . 

Solution. 

By the convolution theorem we have £[h(t)*e~ t *e~ 2t ] = £[/i(f)]£[e -t ]£[e -2 *] = 
s ' ypr ' ypy Using the partial fractions decomposition we can write 

1 1 111 

s(s + l)(s + 2) 2s s + l~^2 s + 2’ 

Hence, 

h(t) * e _t * e~ 2t = ^ — e~ l + ^e _2t , t > 0 ■ 

Problem 47.8 

Compute t * e _t * e l . 

Solution. 

By the convolution theorem we have £[t * e~ l * e*] = £[f]£[e _t ]£[e 4 ] = \ ■ 
U-j- • U_. Using the partial fractions decomposition we can write 

1 _ 1 , 1 1 11 
s 2 (s + l)(s-l) “ “? + 2 ' " 2 ' ITT' 

Hence, 

e * e _t 

t* e t * e* = — t + — - —,t > 0 ■ 


Problem 47.9 


n functions 

, --- A -S 

Suppose it is known that h{t) * h{t) * • • • * h(t) 
stants C and the positive integer n. 


Ct 8 . Determine the con- 


110 



Solution. 


n functions 

We know that £[h(t) * h(t) * • • • * h(t)] = (£[h(t)]) n = -U so that £ _1 [U;] = 
= Ct 8 . It follows that n — 9 and C = ^ ■ 

Problem 47.10 

Use Laplace transform to solve for y(t) : 


sin (t — X)y(X)d\ = t 2 


Solution. 

Note that the given equation reduces to sint * y(t) = t 2 . Taking Laplace 
transform of both sides we hnd = q • This implies Y(s) = = 

| + Jr- Hence, y(t) = £ _1 [ 2 + Jr] = 2 + 1 2 , t > 0 ■ 

Problem 47.11 

Use Laplace transform to solve for y(t) : 

y(t ) - [ e {t ~ x) y(X)dX = t. 

Jo 


Solution. 

Note that the given equation reduces to e* * y(t) = y(t) — t. Taking Laplace 
transform of both sides we hnd Solving for Y(s) we hnd 

y(s) = 8 2(J^ 2 ) • Using partial fractions decomposition we can write 


s — 1 
s 2 (s-2) 


_ i i i 

H + ^ +— 1 —. 

S s 2 (s- 2) 


Hence, 


v(f) - + - + -e 2t , t > 0 


Problem 47.12 

Use Laplace transform to solve for y(t) : 

t*y{t ) = t 2 (l - e _< ). 
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Solution. 

Taking Laplace transform of both sides we find = Jj— p+ip • This implies 
y(s) = | — . Using partial fractions decomposition we can write 


{s + l) 3 


12 1 

s + 1 (s + l) 2 (s + l) 3 


Hence, 


y{t) = 2 - 2 (e~* - 2te~ t 


t 2 

+r ) = 2 


(1 




,t > o ■ 


Problem 47.13 

Solve the following initial value problem. 

y'-y= [ (t-X)e x dX, y( 0) =-1. 
Jo 


Solution. 

Note that y' — y — t * eL Taking Lalplace transform of both sides we find 
sY (-1) - Y = js ■ This implies that T'( s ) = -y 1 ! + sHs-i) 2 ' Usin § 
partial fractions decomposition we can write 

1 _ 2 1 2 1 
s 2 (s-l ) 2 s + s 2 s-l + (s-l) 2 ‘ 


Thus, 

_ 2 1_ _ 2_ 1 _ 2 1_ _ 3 _ 1 

S -l + S + 5 2 S — 1 (s — l) 2 S^S 2 S — 1 (s — l) 2 

Finally, 

y(t) = 2 + t — 3e* + te f , t > 0 ■ 
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Section 48 


Problem 48.1 

Evaluate 

(a) / q 3 (1 + e -i )<5(t — 2)dt. 

(b) f* 2 ( 1 + e _t )<5(t — 2)dt. 

Solution. 

(a) / q 3 ( 1 + e _i )<5(t — 2)dt = 1 + e -2 . 

(b) J^ 2 (l + e -t )<5(t — 2)dt = 0 since 2 lies outside the integration interval ■ 

Problem 48.2 

Let /(£) be a function defined and continuous on 0 < t < oo. Determine 

(f*8)(t) = [ f(t-s)S(s)ds. 

Jo 

Solution. 

Let g(s) — f(t — s). Then 

(f*&)(t)=[ f(t-s)8(s)ds = [ g(s)S(s)ds 
Jo Jo 

=9(0) = f(t) ■ 

Problem 48.3 

Determine a value of the constant to such that J 0 ' sin 2 [7r(t — to)\5(t—\)dt = 

Solution. 

We have 

J sin 2 [ir(t-t 0 )]8(t - \) dt =\ 



Thus, a possible value is when n ((| — t 0 ) = §• Solving for t 0 we find to = 
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Problem 48.4 

If Jj 5 t n 5{t — 2 )dt = 8, what is the exponent nl 

Solution. 

We have /, 5 t n 5{t — 2 )dtt = 2 n = 8. Thus, n — 3 ■ 

Problem 48.5 

Sketch the graph of the function g(t) which is defined by g(t) = j* J Q S 5{u — 
1 )duds, 0 < t < oo. 


Solution. 

Note first that J Q S 8{u — 1 )du = 1 if s > 1 and 0 otherwise. Hence, 


^ ^ { J* h(s — 1 )ds — t — 1, 


if t < 1 
if t > 1 ■ 



Problem 48.6 

The graph of the function g(t) = e at 5(t — t 0 )dt, 0 < t < oo is shown. 
Determine the constants a and t 0 . 


« -► 

2 
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Solution. 

Note that 


It follows that to 

Problem 48.7 

(a) Use the method of integarting factor to solve the initial value problem 

y' -y = h(t), y(°) = o. 

(b) Use the Laplace transform to solve the initial value problem cj>' — (j) — 

m, 0(0) = o. 

(c) Evaluate the convolution (p*h(t ) and compare the resulting function with 
the solution obtained in part (a). 

Solution. 

(a) Using the method of integrating factor we find, for t > 0, 


9(t) = 
2 and a = 


0, 0 < t < t 0 

3 at °, to < t < OO 


-1 


y - y =h(t) 

(e-W =e~ t 
e- l y = - e -t + C 
y = - 1 + Ce l 
y = - l + e*. 

(b) Taking Laplace of both sides we find — $ = 1 or f (s) = U-j-. Thus, 
4>(t) = eL 

(c) We have 

((f>*h)(t)= f e^~ s ^h(s)ds— f e^^ds = — 1 + e* ■ 

Jo Jo 

Problem 48.8 

Solve the initial value problem 

y' + y = 2 + 6(t-l), 2 ,( 0 ) = 0, 0 < f < 6. 

Graph the solution on the indicated interval. 
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Solution. 

Taking Laplace of both sides to obtain sY + Y — 2 + e~ s . Thus, Y(s) = 

i(Ai) + S = f-iTi + S- Hen ^ 
y(t) = 


2 - 2e~*, 

2 + (e — 2)e _t , 


t < 1 
t > 1 l 



Problem 48.9 

Solve the initial value problem 

y" = S(t - 1) - 5(t - 3), 2/(0) = 0, y'( 0) = 0, 0 < t < 6. 

Graph the solution on the indicated interval. 

Solution. 

Taking Laplace of both sides to obtain s 2 Y = e~ s — e _3s . Thus, P(s) = 
^ - s?. Hence, 

y(t) — (t — l)h(t — 1) — (t — 3 )h(t — 3) ■ 
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Problem 48.10 

Solve the initial value problem 

y" - 2y' = 5{t - 1), y(0) = 1, j,'(0) = 0, 0 < t < 2. 

Graph the solution on the indicated interval. 

Solution. 

Taking Laplace transform of both sides and using the initial conditions we 
find 

s 2 Y -s- 2{sY - 1) = e- s . 

Solving for S we find Y(s) = ± ^ \ - £ + £. Hence, 

y(t) = 1 - ^h(t - 1) + ^e 2(t_1) /i(f - 1) ■ 



Problem 48.11 

Solve the initial value problem 

y" + 2y' + y = 5(t~2), y( 0) = 0, y\ 0) = 1, 0 < t < 6. 

Graph the solution on the indicated interval. 

Solution. 

Taking Laplace transform of both sides to obtain s 2 Y — 1 + 2 sY + Y — e~ 2s . 
Solving for P(s) we find F(s) = . Therefore, y(t) = fe _t + (t — 

2 )e-( t_2 )/i(t - 2) ■ 
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Section 49 

Problem 49.1 

Find £[y(£)] where 


Solution. 


y(t) 


d_ 

dt 


e *cos2 1 

0 

t + e* 


£[y0)] = £ 


—e t cos 2t + 2e * sin 2 t 
0 

1 + e* 


s+l 


+ 


(s+l) 2 +4 ^ (s+l) 2 +4 

0 

± + -L- 

s ' s— 1 
3-s 

(s+l) 2 +4 

0 


L 8 + FY J 


Problem 49.2 

Find £[y(i)] where 


Solution. 


y(t) 



£[y(f)] 



Problem 49.3 

Find £ _1 [Y(s)] where 


Y (s) 


s 2 +2s+2 

1 

s 2 +s J 
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Solution. 

We have 


Thus, 


C 


C 

-ir_ 


-Ml 

2 


Ls 2 +2s+2J 


1 


C ~ l 1 
£- 


r_ 




L(s+l) 2 +l 2 J 

Ls s+l-l 


2e 4 sin t 
1 - e -4 


C^[Y{s)} = 


2e 4 sin t 




Problem 49.4 

Find £ _1 [Y(s)] where 



1 

1 

bO 

_1 


r m i 

Y(s) = 

2 0 3 


£[e 24 ] 


1-2 1 


£[sinf] 


Solution. 

We have 



' 1 -1 2 


r m i 

Y(>) = 

2 0 3 


£[e 24 ] 


1-2 1 


£[sinf] 


= C 


1-12 
2 0 3 

1-2 1 



= C 


Thus, 


£ -1 [Y(s)] = 


t 3 — e 24 + 2 sin t 
21 3 + 3 sin t 
t 3 — 2e 24 + sin t 

t 3 — e 24 + 2 sin t 
2t 3 + 3 sin t 
t 3 — 2e 24 + sin t 


Problem 49.5 

Use the Laplace transform to solve the initial value problem 


/ 

y 




o 

o 
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Solution. 

Taking Laplace of both sides and using the initial condition we find 


sY 


5 -4 
5 -4 


Y +[0 ;] 


Solving this matrix equation for Y we find 


Y(>) 


4 


s 2 (s—1) 
s —5 

s 2 (s-l) 


Using partial fractions decomposition we find 

4 _ 4 4 4 

s 2 (s — 1) s 2 + s s- 1 


and 

s — 5 

5 4 4 


s 2 (s - 1) 

s 2 s s — 

Hence, 

Y(a) = 

' _4 _1 4_4_ - 

s 2 s s— 1 

_5 _ i 4_ ±_ 

_ s 2 s s— 1 _ 

Finally, 

y (*) = 

4t + 4 — 4e* 

5t + 4 — 4e f 


Problem 49.6 

Use the Laplace transform to solve the initial value problem 


/ 

y = 


5 -4 
3 -2 


y, y(o) 


3 

2 


Solution. 

Taking Laplace of both sides and using the initial condition we find 


sY- 


3 

2 


5 

3 



Y 
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Solving this matrix equation for Y we find 


Y(>) 


1 

(*-l)(a-2) 


' s + 2 —4 

" 3 ' 

3 s — 5 

2 


1 

(s-l)(s-2) 


3s — 2 
2s - 1 


Using partial fractions decomposition we find 


and 


Hence, 


Finally, 


3s-2 


-1 4 

+ 


(s — l)(s — 2) s — 1 s — 


2s - 1 


-1 3 

+ 


(s — l)(s — 2) s — 1 s 


Y (s') = 


r _i_ _J_ 

s —1 ~ s- 2 

+ -A. 

L s —1 ~ s—2 


y (t) = 


—e t + 4e 2t 
—e* + 3e 2t 


2 


2 ' 


Problem 49.7 

Use the Laplace transform to solve the initial value problem 



1 4 


0 


' 3 " 

/ 

y = 

-1 1 

y + 

3e* 

, y(o) = 

0 


Solution. 

Taking Laplace of both sides and using the initial condition we find 


sY- 


' 3 ' 


1 4 ' 

Y + 

0 

0 

— 

-1 1 

3 

_ s —1 _ 


s - 1 -4 

1 s- 1 


Y = 


3 

3 

s —1 
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Solving this matrix equation for Y we find 


Y(>) 


1 

t-H 

1 

CO 

1_ 


3 

(s—1) 2 +4 

-1 s — 1 


3 

_ s— 1 _ 


(s-l) 2 +4 


3(*-l) + ^ 
0 


Using partial fractions decomposition we find 

3(s-l) + ^ 0 s-1 12 

(s- l) 2 + 4 (s- l) 2 + 4 (s- l)[(s- l) 2 + 4]' 


But 

12 

3 

Hence, 

(s - l)[(s - l) 2 + 4] 

s — 1 

‘ 3 


Y (s) = 

s— 1 
0 

Finally, 


y (t) = 

36* 


0 


Problem 49.8 

Use the Laplace transform to solve the initial value problem 


// 


y 


-3 -2 
4 3 


y, y(0) 



o 

l 


Solution. 

Taking Laplace of both sides and using the initial condition we find 



' 1' 


" 0 ' 


" -3 -2 ' 

s 

0 


1 

— 

4 3 


s 2 + 3 2 

-4 s 2 - 3 
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Solving this matrix equation for Y we find 


Yfs) = 


S 4 —1 


' s 2 - 3 —2 

s 

4 s 2 + 3 _ 

1 


S 4 —1 


s 3 - 3s - 2 
s 2 + 4s + 3 


Using partial fractions decomposition we find 
s 3 - 3s - 2 1 


s 4 - 1 


2^-+ 1 


and 


s 2 + 4s + 3 


s — 1 s 2 + l s 2 + l 


s + 3 


(s — l)(s + l)(s 2 + 1) (s — l)(s 2 + l) s —1 s 2 + 1 
Hence, 

Y(s) = 


s 2 + 1 


_I_ L _i_L_ 

s-l ' s 2 +1 ' s 2 +l 

2 s 1 


Finally, 


y {t) = 


s— 1 s 2 +l s 2 +l 

-e t + 2 cos t + sinf 
2e* — cos t — sin t 


Problem 49.9 

Use the Laplace transform to solve the initial value problem 


' 1 -1 ' 
1 -1 

y + 

' 2 ' 
l 

, y(o) = 

' 0 ' 
1 

, y(0) = 

' o' 
0 


Solution. 

Taking Laplace of both sides and using the initial condition we find 


s 2 Y — s 


0 

1 


Y 


i 

L S 


s 2 - 1 1 

-1 s 2 + 1 


Y = 


2 

s 

l s+ S 
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Solving this matrix equation for Y we find 


Y(>) 



s 2 + 1 

-1 


l 

(Ml co 

s 4 

1 

s 2 -l _ 


J + s J 


y W 



s 5 

J_ 

s 5 



1 + - 
1 ^ 4! 


Problem 49.10 

Use the Laplace transform to solve the initial value problem 


y' = 

' 1 0 o ‘ 
0-11 

y + 

e f 

1 

, y(o) = 

' 0 " 
0 


0 0 2 


-2 1 


0 


Solution. 

Taking Laplace of both sides and using the initial condition we find 


■ 1 

0 

0 " 


i - 

s-l 

0 

-1 

1 

Y + 

1 

0 

0 

2 


S 2 

L J 


s — 1 

0 0 


1 

s-1 

0 

0 

s + 1 -1 

0 s — 2 

Y = 

1 

s 

_2_ 

s 2 - 


Solving this matrix equation for Y we find 


Y(>) 


- i 

0 

0 


1 - 

s— 1 

0 

1 

1 


s-l 

1 

s+1 

(s+l)(s—2) 


s 

0 

0 

1 

s—2 


2 

L s 2 J 


1 


(s-l ) 2 

s(s—2)—2 
s 2 (s+l)(s—2) 
-2 

s 2 (s-2) 
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Using partial fractions decomposition we find 


and 


Hence, 


Finally, 


s ( s ~ 2 ) — 2 1 1 1 1 

s 2 (s + l)(s — 2) s 2 2s 3 s + 1 

-2 1 | 1 1 1 
s 2 (s — 2) “ ^ + 2s ~ 2 ’ s — 2’ 


Y(s) 


+ 


2s 



(S-1)2 
1 _J_ 

3 ' s+1 

_L_1 

2s 2 



s—2 


1 

s—2 


y if) 


te l 

i i i _ Ip-t _ l p 21 
'2 3 6 

t+\-\e 2t 


1 

6 


1 

s — 2 


Problem 49.11 

The Laplace transform was applied to the initial value problem y' 

yi(t) 


y 0 , where y (t) = 


V2 (t) 


, A is a 2 x 2 constant matrix, and yo 


The following transform domain solution was obtained 


Ay, y(0) 
2/i,o 
_ y 2,0 


£[y(t)] = Y(s) 


1 

' s-2 

-1 

1/1,0 

s 2 - 9s + 18 

4 

s — 7 

_ 1/2,0 


(a) what are the eigenvalues of A? 

(b) Find A. 


Solution. 

(a) det(sl — A) = s 2 — 9s + 18 = (s — 3)(s — 6) = 0. Hence, the eigenvalues 
of A are r\ = 3 and r -2 — 6. 

(b) Taking Laplace transform of both sides of the differential equation we 
find 

sY - y 0 = AY 

Y = (si - A)- 1 yo 

Letting s = 0 we find 


Y (0) = -A = — 
v ; 18 


-1 

-7 
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Hence, 


and 


It follows that 


A' 1 = 


1 

18 


2 

-4 


1 

7 


det (A x ) 


1 

18' 


A = (A^ 1 ) -1 


7 -1 
4 2 
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